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imaging  system  and  a  laser  transmitter.  For  the  former,  the  OTF  is 
developed  as  a  measure  of  the  systems  performance,  while  for  the  latter, 
the  antenna  gain  is  analyzed.  Initially,  reference  results  are  developed 
for  a  truly  diffraction  limited  system.  These  results  are  then  extended 
to  the  case  of  an  adaptive  optics  system  in  the  presence  of  turbulence, 
but  with  no  angular  separation  between  the  target  direction  and  the 
reference  beacon  direction.  At  this  level  the  subject  of  adaptive  optics 
correction  for  scintillation  induced  "random  apodization"  is  introduced 
and  is  shown  that  fundamentally  different  random  apodization  corrections 
are  appropriate  for  the  two  types  of  systems.  To  allow  treatment  of  the 
case  where  the  target  and  beacon  directions  are  distinct  and  there  is  a 
potential  anisoplanatism  problem,  previous  propagation  theory  is  extended 
to  'develop  a  new  set  of  statistical  results.  (It  is  found,  inter  alia, 
that  the  log-amplitude:  phase  cross  covariance  has  zero  value) .  Using 
these  propagation  results,  OTF  and  antenna  gain  results  are  developed  for 
adaptive  optics  systems  with  and  without  random  apodization  correction. 

In  examining  these  results  it  is  found  that  there  is  no  fundamental 
difference  (other  than  a  minor  random  apodization  loss  factor)  in  perfor¬ 
mance  whether  or  not  random  apodization  correction  is  provided.  It  is 
found  that  for  both  types  of  systems  the  anisoplanatism  dependence  is 
!  governed  by  the  same  function.  Extensive,  normalized  results  are 
presented  in  graphic  form  for  the  case  of  horizontal  propagation  over  a 
uniform  strength  of  turbulence  path.  Easily  used  approximate  results  are 
developed,  based  in  large  part,  on  the  definition  of  the  isoplanatic  angle 
A  definition  of  is  provided  for  uniform  and  nonuniform  paths. 

In  certain  remote  sensing  applications,  it  is  of  interest  to  determine  the 
power  backscattered  from  air  molecules'  when  a  pulsed  laser  beam  propagates 
through  the  atmosphere.  Our  primary  concern  in  Chapter  2  is  the  evaluation 
of  the  interaction  between  the  laser  pulse  and  the  atmosphere  at  altitudes 
above  the  region  where  turbulence  effects  are  appreciable.  The  analysis 
indicated  that  a  significant  fraction  of  the  laser  pulse  is  backscattered 
into  a  1  m^  detector  if  the  pulse  is  observed  as  it  propagates  from  an 
altitude  of  20  km  to  an  altitude  of  30  km.  Analysis  and  numerical  results 
are  presented. 
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ABSTRACT 


The  work  presented  in  this  report  is  concerned  with  two  distinct 
matters.  Each  of  these  is  covered  in  one  of  the  two  chapters  of  this 
volume.  Chapter  1  treats  the  subject  of  anisoplanatism  in  ordinary 
adaptive  optics.  The  subject  of  Chapter  2  is  the  development  of  an 
estimate  of  pulsed  laser  backscatter  signal  strength  from  the  upper 
atmosphere. 

In  Chapter  1  detailed  analysis  is  presented  for  the  performance 
of  an  imaging  system  and  a  laser  transmitter.  For  the  former,  the 
OTF  is  developed  as  a  measure  of  the  systems  performance,  while 
for  the  latter,  the  antenna  gain  is  analyzed.  Initially,  reference  re¬ 
sults  are  developed  for  a  truly  diffraction  limited  system.  These  re¬ 
sults  are  then  extended  to  the  case  of  an  adaptive  optics  system  in  the 
presence  of  turbulence,  but  with  no  angular  separation  between  the  tar¬ 
get  direction  and  the  reference  beacon  direction.  At  this  level  the  sub¬ 
ject  of  adaptive  optics  correction  for  scintillation  induced  "random 
apodization"  is  introduced  and  it  is  shown  that  fundamentally  different 
random  apodization  corrections  are  appropriate  for  the  two  types  of 
systems.  To  allow  treatment  of  the  case  where  the  target  and  beacon 
directions  are  distinct  and  there  is  a  potential  anisoplanatism  problem, 
previous  propagation  theory  is  extended  to  develop  a  new  set  of  statistical 
results.  (It  is  found,  inter  alia  ,  that  the  log -amplitude:  phase  cross 
covariance  has  zero  value.)  Using  these  propagation  results,  OTF  and 
antenna  gain  results  are  developed  for  adaptive  optics  systems  with 
and  without  random  apodization  correction.  In  examining  these  results 
it  is  found  that  there  is  no  fundamental  difference  (other  than  a  minor 
random  apodization  loss  factor)  in  performance  whether  or  not  random 
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apodization  correction  is  provided.  It  is  fopnd  that  for  both  types  of 
systems  the  anisoplanatism  dependence  is  governed  by  the  same  function. 
Extensive,  normalized  results  are  presented  in  graphic  form  for  the 
case  of  horizontal  propagation  over  a  uniform  strength  of  turbulence  path. 
Easily  used  approximate  results  are  developed,  based  in  large  part,  on 
the  definition  of  the  isoplanatic  angle,  «?0  .  A  definition  of  i>0  is  pro¬ 
vided  for  uniform  and  nonuniform  paths. 

In  certain  remote  sensing  applications,  it  is  of  interest  to  deter¬ 
mine  the  power  backscattered  from  air  molecules  when  a  pulsed  laser 
beam  propagates  through  the  atmosphere.  Our  primary  concern  in 
Chapter  2  is  the  evaluation  of  the  interaction  between  the  laser  pulse 
and  the  atmosphere  at  altitudes  above  the  region  where  turbulence  effects 
are  appreciable.  The  analysis  indicates  that  a  significant  fraction  of  the 
laser  pulse  is  backscattered  into  a  1  ma  detector  if  the  pulse  is  observed 
as  it  propagates  from  an  altitude  of  20  km  to  an  altitude  of  30  km. 
Analysis  and  numerical  results  are  presented. 
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Introduction 
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1. 1 

The  use  of  adaptive  optics  for  correction  of  the  effects  of  atmospheric 
turbulence,  when  controlled  by  the  technique  which  has  become  known  by  the 
name  of  "phase  reversal"  or  "phase  conjugation",  relies  on  the  availability 
of  a  suitable  reference  beacon.1  The  value  of  the  turbulence  induced  per¬ 
turbations  of  the  wavefront  of  this  reference  signal  as  measured  at  the 
aperture  of  the  adaptive  optics  system  is  assumed  to  provide  a  valid  indication 
of  the  turbulence  effects  associated  with  the  propagation  path  of  interest.  If 
the  adpative  optics  is  part  of  a  compensated  imaging  system  then  it  follows 
directly,  while  if  the  adaptive  optics  is  part  of  a  laser  transmitter  system  it 
follows  from  consideration  of  reciprosity2,  that  if  the  adaptive  optics  can  apply 

a  correction  which  corresponds  to  distorting  an  undistorted  wavefront  so  that 

* 

it  is  appropriately  inverse  to  the  reference  signals  wavefront  distortion,  then 

the  adaptive  optics  system  will  perform  in  a  diffraction  limited  manner - 

providing  that  the  propagation  path  for  the  reference  signal  exactly  matches 
the  propagation  path  for  which  diffraction  limited  system  operation  is  desired. 

This  last  point  is  of  some  significance  since  there  are  cases  of  interest 

for  the  application  of  adaptive  optics  systems  for  which  the  angular  separation, 
— * 

d,  between  the  direction  of  the  propagation  path  over  which  the  adaptive  optics 

system  is  to  provide  diffraction  limited  performance,  and  the  direction  of  the 

propagation  path  along  which  the  reference  signal  arrives  at  the  adaptive  optics 

system's  aperture,  is  not  zero.  In  such  cases  we  have  to  ask  whether  or  not 

the  angular  separation,  t?,  is  large  enough  to  be  of  noticeable  consequence  in 

the  sense  of  significantly  degrading  the  performance  of  the  adaptive  optics 

—  > 

system.  If  the  magnitude  of  the  angle,  d,  is  small  enough  that  there  is  no 
significant  degradation,  then  we  say  that  the  angle  "lies  within  the  isoplanatic 


*  The  reader  is  cautioned  here  that  exactly  what  constitutes  the  "appropriate 
inverse"  is  not  uniquely  defined  when  the  power  density  of  the  beacon  signal 
varies  across  the  aperture,  due  to  turbulence  effects.  For  the  compensated 
imaging  system  and  for  the  laser  transmitter,  the  correction  for  such  power 
density  variations  are  distinctly  different,  as  we  shall  see. 


sjc 

patch".  If,  however,  the  angular  separation,  >?,  has  a  large  enough 
magnitude  that  there  is  a  significant  amount  ot  degradation  of  the  per¬ 
formance  of  the  adaptive  optics  system,  then  we  would  say  that  there  is 
an  anisoplanatism  problem. 

In  this  chapter  we  shall  concern  ourselves  with  the  task  of  developing 
a  set  of  theoretical  results  that  will  quantify  these  considerations.  We  shall 

treat  two  classes  of  adaptive  optics  systems - an  imaging  system  and  a 

laser  transmitter.  For  the  former  we  shall  be  interested  in  developing 
a  relationship  between  the  mean  value  of  the  turbulence  degraded  optical 
transfer  function  of  the  imaging  system  and  the  angular  separation,  »>, 
between  the  directions  from  the  center  of  the  system's  entrance  aperture 
to  the  reference  signal's  source  and  to  the  object  to  be  imaged.  For  the  latter 
our  interest  wilt  lie  in  the  task  of  relating  the  turbulence  degraded  antenna 
gain  of  the  laser  transmitter  to  the  angular  separation,  >9,  between  the  di  - 
rection  from  the  center  of  the  system's  exit  aperture  to  the  reference  signal's 
source  and  the  direction  to  the  laser  transmitters'  aim -point. 

A  major  point  of  interest  in  the  analysis  we  shall  be  presenting  relates 
to  the  treatment  of  the  variations  across  the  aperture  of  the  beacon  signal's 
power  density.  In  the  customary  consideration  of  an  adaptive  optics  system 
such  variations  are  ignored.  The  adaptive  optics  applies  a  correction  which 
is  the  inverse  of  the  real  turbulence  induced  phase  shift  sensed  on  the  reference 
signal's  wavefront.  When  d  is  identically  equal  to  zero,  and  if  there  is  no 
power  density  variation  such  operation  of  the  adaptive  optics  will  result  in 


*  The  concept  of  isoplanatism  and  of  an  isoplanatic  patch  was  first  introduced 
in  the  analysis  of  the  resolving  power  of  lenses,  as  a  basis  for  stating  the 
assumption  needed  to  justify  working  with  Fourier  transforms  and  thus 
introducing  OTF  and  MTF  concepts.  The  assumption  was,  in  essence,  that 
the  isoplanatic  patch  size  was  large  enough  that  it  could  be  considered  to  be 
infinite.  The  term  "isoplanatism"  as  it  was  used  in  discussing  a  lens  had 
reference  to  the  size,  shape,  and  orientation  of  the  image  of  a  point  source, 
and  thus  was  only  indirectly  (and  not  uniquely)  related  to  the  wavefront 
aberration.  Nonetheless,  the  terms  "isoplanatism",  "anisoplanatism". 
and  "isoplanatic  patch"  seem  clearly  appropriate  to  denote  the  concepts  we 
have  used  them  to  name  in  this  chapter. 
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diffraction  limited  performance.  However,  if  there  is  some  nontrivial 
power  density  variation  the  performance  will  be  less  than  diffraction 
limited.  The  origin  of  this  degradation  is  perhaps  most  easily  indicated 
by  giving  the  name  "random  apodization"  to  the  power  density  variation. 
Apodization,  i.  e.  ,  the  variation  of  the  the  transmission  across  the  aperture 
of  an  optical  system,  will  effect  the  resolution  of  an  imaging  system  and  the 
antenna  gain  of  a  laser  transmitter.  In  general  the  effect  of  apodization  is  the 
reduction  of  side -lobe  levels  at  the  expense  of  resolution  or  antenna  gain.  When 
the  apodization  is  randomly  "chosen"(by  the  turbulence)  we  may  expect  the 
reduction  in  resolution  or  antenna  gain  to  be  rather  substantial  for  the  amount 
of  apodization,  i.  e.  ,  for  the  amount  of  power  density  variation. 

An  ideal  adaptive  optics  system  should  be  able  to  compensate  for  not 

only  the  real  part  of  the  turbulence  induced  phase  variations  but  also  for  the 

* 

imaginary  part,  i.  e.  ,  for  the  power  density  variations.  When  it  does  this 
the  systems  performance  should,  it  would  seem,  revert  to  the  diffraction 
limited  value,  when  there  is  no  anisoplanatism  problem,  i.  e.  ,  when  equals 
zero.  In  the  following  analysis  we  shall  see  that  this  is  indeed  exactly  the 
case  for  the  adaptive  optics  imaging  system.  For  the  adaptive  optics  laser 
transmitter  we  shall  find  that  performance  cannot,  in  general,  be  made 

exactly  equal  to  the  diffraction  limited  value -  and  that  rather  surprisingly, 

under  some  circumstances  the  performance  can  actually  exceed  the  diffraction 
limited  value  of  antenna  gain.  Moreover,  we  shall  find  that  the  nature  of  the 
random  apodization  correction  is  fundamentally  different  for  the  laser  trans¬ 
mitter  from  what  it  is  for  the  imaging  system. 

In  the  next  section  we  shall  develop  expressions  for  the  diffraction 
limited  performance  of  an  imaging  system  and  of  a  laser  transmitter.  Along 
with  that,  this  section  will  also  consider  adaptive  optics  performance  when 

*  We  consider  it  to  be  outside  the  scope  of  this  work  to  discuss  techniques 
by  means  of  which  adaptive  optics  compensation  for  random  apodization 
could  be  achieved. 


{f  equals  zero,  i.  e.  ,  when  there  is  no  anisoplanatism.  In  doing  this  we 
will  necessarily  also  consider  in  this  section  the  control  rules  required 
for  the  random  apodization  correction  part  of  an  adaptive  optics  system. 

In  the  section  after  that  we  shall  treat  the  statistics  of  propagation  through 
turublence,  developing  formulas  for  the  covariance  of  the  real  and  imaginary 
parts  of  the  complex  phase  variations  induced  by  the  turbulence.  These 
results  will  treat  the  case  when  there  is  a  nonzero  value  for  the  angular 
separation,  $  .  The  section  after  that  will  apply  these  statistical  results 
to  the  task  of  developing  expressions  for  adaptive  optics  imaging  systems 
and  laser  transmitters  with  and  without  random  apodization  correction.  The 
final  sections  will  then  be  concerned  with  the  development  of  numerical  results 
and  conveniently  usabLe  approximations,  and  with  a  discussion  of  these  results. 
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In  this  section  we  shall  consider  the  performance  of  an  ideal  system 
(imaging  or  laser  transmitter)  in  two  senses  of  the  word  ideal.  In  the  first 
sense  we  shall  take  ideal  to  mean  that  there  are  no  turbulence  effects  along 
the  propagation  path,  so  that  system  performance  is  truly  diffraction 
limited.  This  will  be  the  case  treated  in  subsection  1.2.  1  .  Our  objec¬ 
tive  there  will  be  to  present  some  basic  formulations  and  to  develop  some 
reference  results  against  which  other  results  can  be  compared.  Along  the 
way  some  of  the  fundamental  notation  and  assumptions  will  be  introduced. 

In  its  second  sense  of  the  term  "ideal  system"  we  shall  treat  the 
term  as  meaning  that  there  is  no  anisoplanatism  problem,  i.  e.  ,  that  & 
equals  zero.  This  will  be  treated  in  subsection  1.2.2  .  Our  objective 
there  will  be  to  explore  in  a  limited  way  the  implications  of  the  random 
apodization.  It  is  in  this  subsection  that  we  will  develop  an  expression 
for  the  effect  of  random  apodization,  i.  e.  ,  for  beacon  power  density 
variations,  when  there  is  no  provision  in  the  adaptive  optics  for  correct¬ 
ing  for  turbulence  induced  random  apodization.  This  subsection  will  also 
consider  the  question  of  developing  an  appropriate  control  law  for  adaptive 
optics  compensation  of  the  random  apodization.  It  is  in  this  subsection 
that  it  will  be  shown  that  the  control  law  for  random  apodization  correction 
is  fundamentally  different  for  an  imaging  system  than  it  is  for  a  laser  trans¬ 
mitter,  and  that  while  the  imaging  system's  performance  can  be  made  to  be 
exactly  equal  to  the  diffraction  limited  value,  this  is  not  the  case  for  the 
laser  transmitter. 


1.2.1  Performance  in  the  Absence  of  Turbulence  Effects 


In  this  subsection  we  shall  consider  the  performance  of  an  imaging 
system  and  of  a  laser  transmitter  in  the  absence  of  any  turbulence  effects, 
each  assumed  to  be  ideally  fabricated  and  aligned.  For  the  imaging  system 
we  shall  take  the  optical  transfer  function  (OTF)  of  the  system,  as  a  function 
of  spatial  frequency  to  be  our  measure  of  system  performance.  For  the 
laser  transmitter  we  shall  take  the  antenna  gain,  as  our  measure  of  system 
performance. 


Here,  and  throughout  this  chapter  we  shall  assume  that  the  relevant 

wavelength  is  This  will  be  the  same  wavelength  for  the  reference  beacon, 

the  laser  transmitter,  and  the  imaging  system.  Likewise,  throughout  this 

chapter  we  shall  assume  that  our  optical  system  has  a  circular  unobscured 

aperture  of  diameter  D.  Accordingly,  taking  r  to  denote  a  two-dimensional 

position  vector  on  a  plane  containing  the  aperture  with  the  origin  corresponding 

to  the  center  of  the  aperture,  we  can  define  the  aperture  by  the  function 
— ► 

W  (r),  where 


W  (r) 


if  |r|  sD/2 
if  |r|  >D/2 


(2.  1) 


We  shall  consider  a  reference  beacon  located  at  a  range  R  from  the  (center 
of  the)  aperture,  i.  e.  ,  from  the  origin  and  displaced  from  the  z-axis  di¬ 
rection  by  an  angle  t?,  and  a  laser  aim -point  orapoint  source  to  be  imaged 

located  at  the  same  range,  R,  but  displaced  from  the  z-axis  direction  by  an 
,  —* 

angle  +  g  if  .  Thus  the  angular  separation  between  the  point  source  to  be 
imaged  or  the  laser  aim -point  and  the  reference  beacon  is  equal  to  #  .  The 
z-axis  is  nominally  perpendicular  to  the  aperture  plane.  Based  on  the 
assumption  that  |i>|  is  very  small,  then  the  z-axis  defines  the  nominal 
propagation  direction. 

► 

The  path  length  from  a  point  r  in  the  aperture  plane  to  the  beacon 
is  L  (r,  -ftf).  while  that  to  the  point  source  or  aim -point  is  L  (r,  gt>),  where 

L  (r,  0)  =  (R®+  |r-?R|2)1/s  .  (2.2) 
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Based  on  the  assumption  that  |rj  is  much  smaller  than  R,  as  well  as 
that  |  a  |  is  very  small  we  can  rewrite  Eq.  (2)  in  the  approximate  form 

L  (r,  B)  ^  R+  ^  -  r  •  0  +  |  Rfl2  .  (2.  3) 

To  provide  focusing  at  range  R  the  optical  power  elements  of  the  optical 
system  must  produce  an  r  -  dependent  optical  path  length  O  (r)  ,  where 

O  (7)  =  -  |  r2/R  .  (2.  4) 

With  these  notational  matters  taken  care  of  we  are  now  ready  to 
consider  the  development  of  expressions  for  the  diffraction  limited  OTF 
and  antenna  gain.  We  shall  take  up  the  OTF  problem  first. 

1.2. 1.1  OTF  in  the  Absence  of  Turbulence  Effects 

In  this  section  we  shall  fairly  closely  follow  the  imaging  system 
analysis  presented  in  one  of  our  earlier  works3.  We  start  by  noting  that 
if  the  wavefunction  from  a  point  source  [at  (R,  3)  ]  after  passing  through 
the  optics  (and  adaptive  optics,  if  any)  can  be  written  as  the  pupil  function 
U  (r),  then  for  an  optical  system  with  focal  length  F  the  wavefunction, 
u  (x),  at  position  x  in  the  focal  plane  of  our  imaging  system,  where  x  is 
a  two-dimensional  position  vector  on  the  focal  plane,  can  be  written  as 

u  (x)  =  fd?  U  (r)  exp  (-ikr  •  x/F)  ,  (2.  5) 

where  here  and  throughout  this  chapter 

k  =  2  tt/X  •  (2.6) 

The  power  density  in  the  focal  plane  associated  with  this  wavefunction  can 


be  written  as 


9  (x)  =  |  |  u  (x) 


=  v  u  (x)  u  (x)  .  (2.  7) 

Our  interest  is  in  the  OTF,  which  can  be  evaluated  as  the  Fourier 

-4  *4 

transform  of  the  power  density,  9  (x),  since  9  (x)  represents  the  image  of  a 
point  source.  Accordingly  we  can  write  for  the  optical  transfer  function  for 
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spatial  frequency  f  (with  f  measured  in  cycles  per  radian -field -of -view) 

X  (I)  =  B  f  dx  &  (x)  exp  ( -2  n i f  *  x/F)  ,  (2.  8) 

where  B  is  a  normalization  constant  choosen  to  make  the  OTF  of  zero 
spatial  frequency  equal  to  unity,  i.  e.  ,  so  that 

1  (0)  =  1  ,  (normalization  condition)  .  (2.  9) 

When  we  substitute  Eq.  (5)  into  Eq.  (7)  twice,  once  with  the  variable 
of  integration  denoted  by  r  and  once  by  r'*,  and  then  substitute  that  result 
into  Eq.  (8),  we  can  after  making  the  product  of  integrals  into  a  double 
integral,  write  the  result  as 

X  (?)  =  iB  X"2  F“2  Tf  dx  dr  d7'  U*  (r')U  (7) 

X  exp  {  i(k/F)  x  •  [  7'-  (7  +  X?)]  }  .  (2.10) 

We  note  that  the  combination  of  the  x -integration  and  the  r  -  integration  in 
Eq.  (10)  can  be  considered  to  be  a  repeated  Fourier  integral.  As  is  well 
known4,  in  accordance  with  the  equation 

JJ*  dp  *  dg  g(p')  exp  [  ±ig  •  (p'-p)]  =  (2n)n  g(p)  ,  (2.  11) 

where  g(p)  is  any  reasonably  well  behaved  function  of  p  ,  and  p,  p',  and 
g  are  n -dimensional  vectors,  the  repeated  Fourier  integral  recovers  the 
starting  function.  Accordingly  we  can  rewrite  Eq.  (10)  as 

2  (?)  =  £B  X"2  F-2  (F/k)2  (2  tt)2  !  dr  l/  ( r  +  \  f)  U  (r  ) 

=  4B  fd?U*  (r  +  >.f)  U  (r  )  ,  (2.12) 

reducing  our  result  to  a  simple  integral  over  the  aperture. 

To  evaluate  the  normalization  constant,  B  ,  we  make  use  of 
Eq's.  (9),  and  (12).  This  allow  us  to  write 

B  =  {  TdT  [|  U*(7)  U  (r)]  I-1  •  (2.  13) 

•J 
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r 


Further  simplification  of  this  result  is  generally  possible  based  on  the 
fact  that  in  this  case  the  power  density  in  the  aperture,  |  U  (r)U(r), 

and  in  the  more  general  case  the  average  power  density  in  the  aperture, 

— ♦  — *  — , 

<  2  U  (r)U(r)),  is  not  a  function  of  aperture  position,  r  .  Thus  we 

can  write  for  the  power  density  in  the  aperture 

=  |u*(?)  u(7)  ,  (2. 14) 

or  in  the  more  general  case 

e,  =  <  |U*(r)  U(?)  y  ,  (2.  14') 

where  the  angle  brackets,  (  .  .  .  )  ,  denote  an  ensemble  average.  Making 
use  of  Eq.  's  (1)  and  (14),  Eq.  (13)  can  be  reduced  to  the  form 

B  =  J*  dr  W(r)  r1 

=  (i  tt  D2  )-1  .  (2.15) 

For  the  diffraction  limited  case  which  we  are  concerned  with  here 
— ► 

the  pupil  function  U  (  r  )  is  defined  by  an  amplitude,  A  ,  which  is  not  a 

— ► 

function  of  position,  r  ,  a  path  length  induced  phase  shift,  L  ,  and  a 
focusing  phase  shift, O.  Thus  we  can  write  for  a  source  at  (  R,  6  ) 

U  (  r  )  =  A  W  (  r  )  exp  fik[L(r,0)  +  O(r)]],  (2.  16) 

where  we  recall  that  L,  O,  and  k  are  as  defined  by  Eq.  's  (2)  or  (3),  (4), 
and  (6),  respectively.  We  can  rewrite  this  as 

U(7)  =  A  W  (7)  exp  [  ik(R  +  #  R  02  -r-0)]  ,  (2.17) 

from  which  in  conjunction  with  Eq.  (14),  it  follows  that 

9k  =  i  A*  A  W  (7)  ,  (2.  18) 

and 

B  =  (-^  ttD2  A*  A  )_1  .  (2.  19) 

When  we  substitute  Eq.  's  (17)  and  (19)  into  Eq.  (12)  and  carry  out 


the  obvious  simplifications  we  get  the  result  that 


(2.  20) 


r 


2(f)  =  (iuD2)-1  Jdr  W(r)W(rUf  )  exp(ik  \7  •  6  ) 

=»  ( i  tt  D2  T1  Jdr  W  (  r  )  W  (  r.  +  X? ) 

— *  — ► 

In  writing  Eq.  (20)  we  have  dropped  the  factor  exp  (ikxf  •  8)  as  this  is 
just  a  phase  shift  factor  indicative  of  the  fact  that  the  point-source  is  not 
on  the  z -axis  but  displaced  from  it  by  an  angle  G.  The  modulus  of  this 
factor  is  unity  and  for  the  on -axis  point-source  the  factor  would  be  exactly 
unity.  Thus  in  dropping  this  factor  we,  in  essence,  make  the  point -source 
define  the  on -axis  direction  for  OTF  purposes. 

The  integral  in  Eq.  (20)  represents  the  area  of  overlap  of  two  circles 
of  diameter  D  with  centers  displaced  by  amount  K  f .  It  is  easy  to  show 
from  some  simple  trigonometric  analysis  that  if  two  circles  of  diameter  D 
have  their  centers  separated  by  a  distance  p  ,  then  the  area  of  overlap  of 
the  circles  is  4ttD2K(p),  where 

.f-  [  cos-1  (o/D)  -  (p/D)  [l-(p/D)s]1/a  J  ,  if  p<  D 
K(p)  =] 

U  ,  if  p  >  D  .  (2.  21) 

Thus  we  can  rewrite  Eq.  (20)  as 
Xol  (?)  =  K(\f) 

%  cos-1  (Xf/D)  -  (Xf/D)  [1  -(Xf/D)s]1/s  }  ,  if  XfsD 
{o  ,  if  Xf  >  D  .  (2.  22) 

The  subscript  "DL"  is  added  here  to  make  explicit  the  fact  that  this  result 
applies  for  the  case  of  diffraction  limited  operation.  Eq.  (22)  represents 
our  basic  result  for  the  OTF  of  an  ideal  imaging  system  in  the  absence  of 
turbulence.  In  the  next  subsection  we  shall  consider  the  corresponding 
problem  for  a  laser  transmitter. 

1.2.  1.2  Antenna  Gain  in  the  Absence  of  Turbulence 

The  antenna  gain  of  a  laser  transmitter,  G  ,  is  most  usefully 
defined  as  the  ratio  of  the  laser  powered  density  at  the  aim -point,  , 

divided  by  the  total  laser  power  transmitted,  PT  ,  and  scaled  for  the  range, 
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R  ,  from  the  laser  transmitter  to  the  aim -point.  Thus  we  would  write 


G  =  (0AP  / PT  )  R2  .  (2.  23) 

Expressed  in  this  form  the  antenna  gain  has  the  dimension  of  inverse 
steradians.  It  is  in  effect  a  measure  of  the  (inverse  of  the)  beam  spread 
of  the  laser  transmitter. 

If  the  laser  wavefunction  leaving  the  aperture  is  U  (r),  then  the 

wavefunction  at  an  aim -point  at  (R,  9  )  can  be  written  as 

uAP  =  -  -  Jdr  U  (r)  exp  [  ikL  (  r  9  )  ],  ,  (2.  24) 

where  L  is,  as  defined  by  Eq.  's  (2)  or  (3),  the  path  length  between 
— ♦ 

the  point  at  r  on  the  aperture  plane  and  the  aim -point.  To  focus  the  laser 
beam  on  the  aim -point  the  laser  transmitter  optics  would  cause  the  wave - 
function  leaving  the  aperture  to  have  a  form  expressable  as 

U(7)  =  A  W(r)  exp  {ik[0(O  +7-  9  ]  }  ,  (2.25) 

— ►  — * 

where  the  r  •  U  -  term  in  the  exponential -function  indicates  that  the  laser 
beam  is  directed  at  the  aim -point,  and  the  O  (r)-term  indicates  that  the 
laser  beam  is  focused  at  the  range,  R  ,  of  the  aim -point.  Making  use  of 
Eq.  's  (3),  (4),  and  (25),  Eq.  (24)  can  be  reduced  to  the  form 

uAP  =  exp  [ikR(l  +  i92)]  A  J*  d7  W  (r)  .  (2.26) 

Taking  note  of  Eq.  (1)  we  can  further  reduce  this  to 

Uap  =  -y-R  exp  [ikR(l+|93)]  A  |ttDs  .  (2.27) 

The  laser  power  density  at  the  aim -point  can  be  written  as 

i  * 

^»p  -  5  uir  u»f  .  (2.  28) 

*  The  antenna  gain  as  we  have  defined  it  can  be  related  to  the  more  con¬ 
ventional  antenna -gain -relative -to -isotropic,  by  multiplying  our  value 
by  4  rr. 


(2.  29) 


which  can  be  reduced,  by  means  of  Eq.  (27)  to  the  form 

=  k  R~3  A*  A  &nDs)? 

The  Laser  power  density  at  the  aperture  is  as  defined  by  Eq.  (14)  and 
in  view  of  Eq.  's  (1)  and  (2  5)  can  be  written  as 

<?A  -  i  A*  A  W  ("r )  .  (2.  30) 

The  total  laser  power  transmitted  is  of  course  just 

PT  =  j  dr  9 A 

-  |  A*  A  J‘dr  W  (  r  ) 

=  ttD2  A*  A  .  (2.  31) 

When  we  substitute  Eq.  's  (29)  and  (31)  into  Eq.  (23)  we  obtain  for 
the  antenna  gain,  the  result  that 

Gol  =  in(D/\)2  .  (2.32) 

Here  again  the  subscript  "DE"  has  been  added  to  make  explicit  the  fact 
that  this  result  applies  for  the  case  of  diffraction  limited  operation. 

With  these  results  in  hand,  i.  e.  ,  Eq.  's  (22)  and  (32)  we  have  com¬ 
pleted  the  desired  analysis  of  system  performance  in  the  absence  of  tur¬ 
bulence.  In  the  next  subsection  we  shall  turn  our  attention  to  the  other 
sense  in  which  we  used  the  term  "ideal  system",  namely  an  adaptive  optics 
system  in  which  there  is  no  angular  separation  between  the  directions  from 
the  transmitter  to  the  point -source  or  aim -point  and  to  the  reference  beacon. 


1.  2.2  Adaptive  Optics  Performance  With  Turbulence  But  No  Angular 


Separation 


To  make  our  notation  sufficiently  general,  so  that  it  can  be  used  in  the 

latter  parts  of  this  chapter,  we  shall  consider  a  reference  beacon  located  at 

— * 

(R,  08)  and  a  target  (either  a  point-source  to  be  imaged  or  a  laser  aim -point) 
— * 

located  at  (R,  0T  ) .  In  this  section  we  will,  at  an  appropriate  point  introduce 

-4 

the  fact  that  0B  =  0  and  that  0T  =  0  ,  but  for  much  of  the  deviation  we  shall 
ignore  these  facts  so  as  to  keep  the  results  generally  applicable.  We  shall 
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consider  an  adaptive  optics  system  which  can  perfectly  sense  the  wave  - 
function  perturbation  induced  on  the  beacon  signal  by  propagation  through 
turbulence,  and  which  can  apply  exactly  the  desired  correction.  With  no 
angular  separation  between  the  directions  there  is  no  reason  why  the  system 

performance  should  not  be  ideal  -  except  for  the  matter  of  how  to  handle 

power  density  variations  across  the  system's  aperture.  And  it  is  indeed 
just  this  matter  that  is  the  central  concern  of  this  subsection. 


Before  we  start  a  detailed  examination  of  this  matter  it  is  necessary 

that  at  this  point  we  introduce  the  appropriate  notation  for  the  propagation 

statistics.  For  a  point  source  located  at  position  r  on  the  aperture  and 

traveling  to  an  "end-point"  position  defined  by  the  range  R  and  angular 

deviation  from  the  z -axis  of  0,  i.  e.  ,  by  (R,  0),  if  the  wavefunction  at  the 

end-point  should  be  Uq  (  r,  0)  if  there  were  no  turbulence  along  the  propagation 

— ^ 

path,  and  if  the  instantaneous  random  value  of  the  wavefunction  is  u(r,  0), 
it  is  convenient  to  make  manifest  the  turbulence  effects  by  writing 


u  (  r, 0)  =  u0  (  r,  0)  exp  [i t  (  r, 0)]  (2.  33) 

— >  — > 

The  quantity  h  (  r,  0),  which  we  call  the  "complex -phase  perturbation"  serves 

— >  — > 

as  a  measure  of  the  effects  of  turbulence.  It  is  convenient  to  write  i|i  (  r,  0)  in 

— *  — *  — * 

terms  of  its  real  and  imaginary  parts,  0  (  r,  0)  and  -£(r,0)  respectively.  Thus 
we  have 

*  (r,0)  =  0  (r,0)  -  £  (r,0)  .  (2.34) 

The  quantity  0  (  r,  0  )  is  called  the  "phase  perturbation"  or  the  "real  phase 

perturbation".  It  is  generally  what  is  intended  when  one  speaks  of  "correcting 

—  ►  ► 

the  phase"  of  a  turbulence  distorted  wavefunction.  The  quantity  i(r,  0  )  is 
generally  called  the  "log -amplitude  perturbation".  It  is  a  (logarithmic) 
measure  of  power  density  variations. 

By  virtue  of  the  reciprocity  theorem2  the  complex  phase  perturbation  i 

ijt  (  r,  0  ),  is  not  only  a  measure  of  the  turbulence  induced  perturbation  of  a  point 

source  traveling  from  position  r  on  the  aperture  plane  to  a  target  position  at 

(R,  0),  but  serves  equally  well  as  a  measure  of  the  turbulence  induced  point  source 

-*  ^  _  — * 

at  the  target  position  (R,0)  propagating  back  to  the  position  r  on  the  aperture 
plane.  As  a  consequence,  l  (  r,  0e)  is  a  measure  of  the  power  density  variation  of  the 
beacon  signal  at  the  aperture. 
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The  key  to  our  analysis  in  this  subsection  lies  in  our  ability  to 
form  the  ensemble  average  of  the  exponential  function  of  some  linear 
combination  of  complex  phase  terms.  The  key  to  being  able  to  do  this 
lies  in  the  fact  that  if  x  is  a  gaussian  random  variable,  or  linear 
combination  of  several  gaussian  random  variables,  then  it  is  easy  to 
show  by  carrying  out  an  integration  over  the  probability  density,  that 

<  exp  (ay))  =  exp  (a  X)  exp  [§  a8  (  (  x  "  X  f  )  1  >  (2.  35) 

where  a  is  a  constant,  and 

~X  =  (  X  )  »  (2-  36) 

is  the  mean  value  of  our  random  variable(s). 

Use  of  Eq.  (35)  allows  proof,  by  means  of  conservation  of  energy 
arguments  ,  that  for  the  log  -amplitude  pertubation,  JL  , 

I  =  <  (  i  -  If  )  .  (2-  37) 

— ►  — ♦ 

where  we  have  written  i  in  place  of  l(r,  0)  simply  as  a  matter  of 

convenience.  With  these  results  in  hand  we  are  now  ready  to  take  up  the 

first  of  the  several  cases  to  be  considered  in  this  subsection,  namely  the 

ensemble  average  OTF  of  an  adaptive  optics  imaging  system  when  there  is 

— > 

no  separation  between  the  point -source  (target)  direction,  0T  ,  and  the 

reference  beacon  direction,  0B  ,  for  the  case  of  adaptive  optics  that  only 

1  ¥ 

corrects  for  the  real  phase  perturbations,  0  (  r,  0),  and  not  for  the  power 
density  or  log -amplitude  perturbations,  JL  (  r,  0),i.  e.  ,  not  for  random  apodization. 

1.2.  2.1  Ideal  Adaptive  Optics  Imaging  But  With  No  Random  Apodization 
Correction 

Our  starting  point  for  the  OTF  analysis  to  be  performed  here  con¬ 
sists  of  Eq.  's  (12)  and  (13)  in  ensemble  average  form.  The  ensemble  aver¬ 
age  OTF  (which  can  be  shown  to  correspond  both  to  the  ensemble  average  of 
the  instantaneous  random  OTF,  and  to  the  OTF  associated  with  an  ensemble 
average  of  the  instantaneous  random  image)  can  be  written  as 
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(2.  38) 


<  T(7)  >  =  i  (B  >  fdr  ('u*(rU£)U(?)) 
with  the  ensemble  average  normalization  constant  written  as 

<  B  '  =  f  fd7  £  (  U*  (7)  U(7)>  }-x  ,  (2.39) 

It  is  to  be  noted  that  <  B  )  is  not  the  ensemble  average  of  the  instantaneous 
random  values  of  the  normalization  constant  B  ,  and  thus  is  not  obtained 
directly  from  Eq.  (13)  by  applying  the  ensemble  average  brackets,  ) 
to  both  sides  of  Eq.  (13).  Rather  the  value  of  (  B  \  is  obtained  from 
Eq.  (38)  in  just  the  same  way  that  Eq.  (13)  is  obtained  from  Eq.  (12),  i.  e.  , 
by  reference  to  Eq.  (9),  the  relevant  form  at  which  we  write  here  as 

<  Z  (0)  >  =  1  .  (2.  40) 

Eq.  's  (38)  and  (30)  would  be  relevant  to  any  imaging  system - 

its  just  a  matter  of  the  nature  of  the  specifications  of  the  pupil  function, 

— 

U  (  r  )  .  To  make  it  explicit  here  that  we  are  dealing  with  an  adaptive  optics 
imaging  system  without  any  random  apodization  correction  we  shall  use  the 
subscript  IW/O.  Thus  we  would  write  as  the  locally  relevant  forms  of 
Eq.  's  (38)  and  (39) 

<  'Xlw/o(n  >  =  i  <  B,  w/0  >  f  d"r  .<  U*lw/0  (7+\7)Ulw/0  (7)  >  ,  (2.  41) 

^  m  / o  ^  ~  (  :  h  s  (  U  m/j  (r  )U|i/o  (i)  )  ]  .  (2.  42) 

The  pupil  function,  Uj  w/  0  (  r  )  >  is  in  part  as  given  by  Eq.  (16)  ,  but  in 

addition  there  is  an  exponential  function  of  the  complex  phase,  i]i  (  r,  0T  )  , 

associated  with  propagation  through  turbulence  from  the  (target)  point-source 
— ^  ^  ^ 

to  be  imaged  at  0T)  less  the  turbulence  induced  real  phase  0  (  r,  S8  ),  sensed  in 

viewing  the  reference  beacon  at  0B  .  Thus  we  can  write 

U|  w/o  (  r  )  =  A  W  (  r  )  exp  fi  k  [L  (  r,  0T  )  +  G(  r  )  ] 

+  i  [  4  (  r,  0T  )  -  0(  r,  0B  )  ]  }  .  (2.  43) 
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Making  use  of  Eq.  's  (3),  (4),  and  (34)  this  can  be  rewritten  as 

U,  h/  o  (  r  )  =  A  W  (  r  )  exp  {ik  [R  (1+|9T  )  +  0(  r  )  1 

+  i  [*(?,  eT)  -0(r,  0B)]  }  ■  (2.  44) 

When  we  twice  substitute  this  expression  into  Eq.  (42)  and  appropriately 
simplify  the  results  we  get 

<  BIW/0)  =  a"a  rd7  W(7)(  exp  [2X(7,  0T)  ]  >  f1  .  (2.45) 

Stationarity  of  the  statistics  of  the  real  and  imaginary  parts  of  the  complex 
phase  insures  that  neither  the  mean  value  of  the  log -amplitude 

l  =  it  (7,  9)  >  ,  (2.  46) 

nor  the  covariance  of  the  log -amplitude 

C*(7,£)  =  <[X(7+7.  0  +“3)  -I  (7,  6)  -  I]  )  ,  (2.47) 

— ►  -4 

is  a  function  r  or  0  .  In  accordance  with  Eq.  (35)  we  can  thus  rewrite 
Eq.  (45)  as 

iBIM/0>  =  [|  A*  A  exp  (2  I)  exp  [  2  C£  (0,  0)]  J  dr  W  (  r  )  )-\  (2.48) 

By  virtue  of  Eq.  (37)  the  log -amplitude  mean  value,  £,  and  variance, 
C*(0,0),  exactly  cancel.  [This  exact  cancellation  is  directly  associated 

with  the  fact  of  energy  conservation  -  which  is  what  gave  rise  to  Eq.  (37) 

in  the  first  place.  ]  Thus  we  can  write 

(  Bin/0  >  =  U  A*  A  '  dT  W  (?)  ]"x  ,  (2.  49( 

which  in  view  of  Eq.  (1)  can  be  reduced  to  the  form 

<  B,  1  =  (|  rrD2  A*  A  )“l  .  (2.50) 

This  result,  which  is  identical  to  Eq.  (19)  for  the  no  turbulence  case, 
applies  whether  or  not  there  is  any  angular  separation  between  the  target 

— 4  «4 

point -source  direction,  9T  ,  and  the  reference  beacon  direction,  3e  .  For 

the  balance  of  this  section  we  shall,  however,  restrict  our  attention  to  the 

more  limited  (ideal)  case  which  is  the  subject  of  this  subsection,  where 
— ♦  — * 

9t  and  9B  are  both  equal  to  zero,  so  that  both  target  point -source  and 
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reference  beacon  lie  on  the  z  -axis. 


With  this  restriction  Eq.  (44)  reduces  to 

UIW/Q  (  r  )  =  A  W  (  r  )  exp  [i  k  R  +  X  (  r ,  0)  ]  .  (2.  51 ) 

From  this  it  follows  that 

<  U*lw/0  (r  +  7)  UtV0  (7)  >  =  A*  AW  (7  +  7)  W(7) 

X  <  exp  [2  (  r  +  p,  0)  +  2  (  r,  0  )]>  .  (2.52) 

Making  use  of  Eq.  (52)  to  allow  evaluation  of  the  ensemble  average  of  an 

exponential  function,  we  can  easily  show  that 

(  exp  [ X  (  r  +  p,  0)  +  2  (  r,  0)  ]  )  =  exp  (2  2) 

XexpCCx(0,0)  +  Ci  (7,  0)  ]  .  (2.53) 

F rom  Eq.  (37)  it  then  follows  that 

(exp  [2  (7  +  7,0)  +  2(7,0)])  =  exp[-Cx(0,0)  +  Cx  (7,0)].  (2.54) 
When  we  substitute  this  expression  into  Eq.  (53)  we  get  the  result  that 
(  U*,  w/  o  (  r  +  p  )  U,  w/  o  (  r)  >  =  A*  A  W  (  r  +  p  )  W  (  r  ) 

X  exp  [-C^O,  0)  +  C£(7,  0)]  .  (2.  55) 

By  means  of  this  result  we  can  now  reduce  Eq.  (41)  to  the  form 

S|*/o(7>  =  |A*A  exp'E-Cj  (0,  0)  +  CjUf^O)  ]  <BIH,o  > 

xfdrW(r  +  Xf)W(r)  .  (2.  56) 

This  integral  over  r  is  the  one  discussed  just  before  Eq-  (22).  Making 
use  of  this  fact  we  can  write 

T,w/o(?)  =  |  A*  A  expC-Cx(0,0)  +  Ci(\7,0)]  (  BtH/0  ) 

X  \  ttD3  K  (  Xf )  ,  (2.  57) 

which,  in  view  of  Eq.  's  (22)  and  (50)  reduces  to 
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(2.  58) 


3-ih/o  (?)  =  exp  [  -C^  (0,  0)  +  C_j  (  X?,  0)  ]  K  (  Xf  ) 

=  exp  C  -Cx  (0,  0)  +  Cx(Xf^0)]  X0L  (?) 

The  key  thing  to  note  in  interpreting  this  result  is  that  for  very  low  spatial 
— ♦ 

frequencies  C %  (  Xf ,  0)  is  very  nearly  equal  to  (0,0),  so  that  the  adaptive 

optics  imaging  system  without  random  apodization  correction  provides 

— 

essentially  the  diffraction  limited  OTF,  XDL  (f  )  - while  for  the  higher 

spatial  frequencies,  which  are  what  are  really  of  interest  to  us,  CU  (  Xf  ,  0) 
is  very  nearly  equal  to  zero,  so  that  the  OTF  is  less  than  diffraction  limited 
by  a  factor  of  exp  [  -C  ^  (0,  0  )  ]  .  If  the  log  -amplitude  variance  is  small 
enough  this  factor  is  very  close  to  unity  and  there  would  seem  to  be  no  real 
need  for  a  random  apodization  correction  capability  in  the  imaging  system's 
adaptive  optics.  If,  however,  the  log -amplitude  variance,  C ^  (0,15),  is  not 
particularly  small  then  the  exponential  factor  can  result  in  a  nontrivial  re¬ 
duction  factor - so  that  it  might  be  useful  to  be  able  to  have  the  adaptive 

optics  compensate  for  random  apodization.  We  take  up  an  analysis  of  this 
case  next. 

1.2.  2.  2  Ideal  Adaptive  Optics  Imaging  With  Random  Apodization  Correction 

Using  the  subscript  IW  to  denote  the  fact  that  the  expression  relates 
to  an  adaptive  optics  imaging  system  with  compensation  of  random  apodization, 
we  would  write  in  place  of  Eq.  's  (41)  and  (42),  the  equations 

(  3*1  h  (?)  >  =  i<  B,„  >J  dr  (  U*  H  (7  +  X?)UIH  (7)  >  ,  (2.  5  9) 

for  the  ensemble  average  OTF,  and 

<B,„>  =  {  [dr  *<  U*  H  (7)Uim  (7))  r1  ,  (2.60) 

for  the  cor  responding  normalization  constant. 

For  adaptive  optics  imaging  when  there  is  no  angular  separation 
between  the  target  point -source  to  be  imaged  and  the  reference  beacon,  it 
is  obvious  that  there  will  be  no  residual  turbulence  effects  if  the  adaptive 
optics  random  apodization  correction  makes  the  system's  transmission 
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inversely  proportional  to  the  power  density  of  the  beacon  reference  signal 
over  the  aperture  .  We  shall  shortly  sec  from  our  equations  that  this  is 
exactly  the  case.  We  start  by  noting  that  this  random  apodization  scheme 
would  come  about  if  the  adaptive  optics  imposed  a  complex  phase  correction 
in  which  the  imaginary  part  were  the  negative  of  the  imaginary  part  of  the 
complex  phase  associated  with  the  beacon.  Thus,  in  place  of  Eq.  (43)  we 
would  have 

U,  w  (  r  )  =  A  W  (  r  )  exp  fik  [L  (  r,  eT  )  +  0(  r  )  ] 

+  i[  v  (r,0T)  J  -  i(r,?B)  ]  ,  (2.61) 

and  in  view  of  Eq.  (34),  and  using  Eq.  's  (3)  and  (4)  this  can  be  rewritten 
as 

U(W  (T)  =  AW  (r  )  exp  |ik[R(l  +  ieTS)  -7  ■  0T  ] 

+  i[<J(7,0T)  -  0  (7,0B)]  +  Utr."^)  -x(r,3B)]}.  (2.62) 

For  the  ideal  case  of  interest  to  us  here,  in  which  there  is  no  angular 
separation  between  the  target  point -source  to  be  imaged  and  the  beacon 
reference,  with  both  lying  on  the  z-axis  so  that  0T  equals  0B  ,  then  obviously 
the  two  complex  phase  function  in  Eq.  (62)  will  cancel  and  we  get 

U,„(7)  =  AW(7)  exp  fik[R(l  +  t0Ts)  -7'8t]  }  .  (2.63) 

This  expression  is  manifestly  independent  of  any  turbulence  effects,  and 
by  comparing  of  it  with  Eq.  (16)  we  can  see  that  it  will  lead  to  the  con¬ 
clusion  that  the  OTF  will  be  diffraction  limited.  If  we  carried  out  the  same 


#  It  is  not  entirely  clear  exactly  how  the  random  apodization  correction 
would  be  affected.  Transmission  reduction  is  no  special  problem,  but 
where  the  beacon  power  density  is  very  large  the  transmission  would  have 
to  be  made  correspondingly  large.  Leaving  out  the  concept  of  laser 
amplifiers  (and  nearly  monochromatic  imaging)  it  would  seem  that  affecting 
random  apodization  correction  would  require  acceptance  of  a  very  low 
average  transmission,  so  that  increases  as  well  as  decreases  in  trans¬ 
mission  would  be  possible  and  we  would  have  to  accept  a  limited  dynamic 
range  for  the  compensation.  In  the  analysis  we  have  considered  this  limited 
range,  as  well  as  the  reduced  average  transmission  as  things  we  could  ignore. 
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calculations  as  in  subsection  1.2.  1.  1  we  would  get  exactly  the  same 
results,  which  we  can  write  as 


(!,-(?))  =  Xdl  (?)  -  (2.64) 

This  very  nice  result  would,  of  course,  not  be  expected  to  hold  if  there 
were  a  nonzero  angular  separation,  >V  ,  between  the  target  direction  and 

the  beacon  direction -  but  this  is  a  matter  that  we  shall  take  up  later.  In 

the  next  two  subsections  we  shall  turn  our  attention  to  the  adaptive  optics 
laser  transmitter  problem  in  the  presence  of  turbulence  effects,  but  for 
the  ideal  conditions  case  of  no  angular  separation  between  the  directions 
to  the  target  aim -point  and  to  the  reference  beacon. 


In  exactly  the  same  way  that  we  arrived  at  Eq.  (23)  as  the  ap¬ 
propriate  expression  for  the  non -random  (or  instantaneous)  laser  trans¬ 
mitter  antenna  gain,  we  would  write  for  the  ensemble  average  antenna  gain 

<G>  =  (  <  )  /  (  P  >  )  Ra  ,  (2.  65) 

where  (  9  )  denotes  the  ensemble  average  power  density  at  the  aim  -point 
and<  P  )  denotes  the  ensemble  average  laser  transmitter  power.  It  is  sig¬ 
nificant  to  note  that  the  ensemble  average  antenna  gain,  (  G  ),  is  in  general 
not  the  ensemble  average  of  the  instantaneous  random  antenna  gain.  i.  e.  ,  Eq. 
(65)  is  not  obtained  by  simply  applying  ensemble  averaging  to  both  sides  of 
Eq.  (23)  .  Rather,  Eq.  (65)  is  presented  on  the  same  basis  as  was  Eq.  (23),  only 
this  time  as  the  ratio  of  average  power  density  and  average  transmitted  power. 

The  ensemble  average  laser  power  density  at  the  aim -point,  <  Q  , 

can  be  written  as  Am 

<  9  >  =  }  <  u  '  u  >  ’  ' 

which  can  be  directly  inferred  from  Eq.  (28).  Here  u  is  the  laser  trans¬ 
mitter  produced  random  wavefunction  at  the  aim -point.  This  wavcfunction 
can  be  written  in  terms  of  the  laser  transmitter’s  random  pupil  function  since 
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it  incorporates  the  effects  of  the  adaptive  optics  corrections,  which  are 
random  in  the  sense  that  they  follow  turbulence  effects.  In  place  of 
Eq.  (24)  we  now  have  the  expression 

u  :  l'drU(r)  exp  {i  [k  L  (  r,  9T  )  4  ,4  (  r,  0T  )  ]  }  .  (2.  67) 

Thus  far  over  notation  has  been  quite  general,  not  indicating  whether 
we  were  allowing  for  random  apodization  correction.  To  write  down  an  ex¬ 
pression  for  the  pupil  function  we  must  be  explicit  on  that  matter.  Using  the 
subscript  LT  W  /O  to  denote  an  adaptive  optics  laser  transmitter  without 
any  random  apodization  correction  we  would  write,  as  an  extension  of  Eq.  (25) 

f'lrH/o  (  r  )  -  A  W  (r  )  exp  fik  [  0(  r  )  +  r  •  9r  ]  -  i0  (  r,  9e  )  ),  (2.  68) 

as  the  adaptive  optics  corrected  pupil  function.  With  the  same  subscript 
notation  Eq.  's  (65),  (66),  and  (67)  would  be  rewritten  as 

(  G(.  rw/o  '*  ”  ((  M  n/i)  )  /  (PtTw/o  1  >  (2.69) 

(  ^un/o  '  1  /  u  lth/o  uLTm/o)  >  (2.  70) 

and 

ut  r  h/o  ~  -  ^  fdr  Uuk/o  (r)  exp  fi[k  L  (  r ,  UT  )  +■  *  (  r ,  9t  )  ]  1.  (2.71) 

Making  use  of  Eq.  's  (3),  (4),  (34),  (67>,  and  (68)  we  can  rewrite  Eq.  (71)  as 
u,.T  w/o  -  -  ^  [i  k  R  (i  +  NrS)]  A  f  dr  W  (  r  ) 

x  exp  f  i[$  (  r,  9t  )  -C(r,O8)]+J(r,0Tp  .  (2.  72) 

Wr  note  that  the  ensemble  average  transmitter  laser  power  can  be  written  as 
/6>irw/'i  '  ^dr^{8iTK/o(r)  U 1  tm/  0  (  r )  1  .  (2-73) 

and  that  from  Eq.  (68)  it  then  follows  that 

<  p. TW/O  >  =  r  ci7  h  A*  A  W  (7) 

=  Ti  D3  A*  A  .  (2.  74) 

The  results  up  to  this  point  have  been  developed  without  making  any  use  of 
the  fac  t  that  the  target  aim -point  and  the  reference  beacon  both  lie  on  the 
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z-axis,  so  that  9r  and  3B  are  both  equal  to  zero  and  there  is  no  angular 
separation  between  the  directions  to  the  target  aim -point  and  to  the  reference 
beacon. 

■  ¥ 

When  we  make  use  of  the  zero  values  of  0T  and  0e  ,  Eq.  (72)  reduces 
to 

Ultimo  =  exp(ikR)Ajdr  W  (  r  )  exp  [  X(  r ,  0)  ]  .  (2.75) 

"We  now  substitute  Eq.  (75)  twice  into  Eq.  (70),  make  a  double  integral  out 
of  the  product  of  integrals,  and  interchange  the  order  of  integrations  and 
ensemble  averaging.  Thus  we  obtain  the  result  that 

<  )  =  iX*  R'2  A*  //dr  d?'W(r)  W  (7') 

x(exp  [X(r,  0)  +  X(7',0)  .  (2.76) 

It  is  convenient  at  this  point  to  introduce  the  variable 

p  =  t*  -  r  ,  (2.  77) 

which  on  making  a  change  of  variable  of  integration  in  Eq.  (76)  allows  us 
to  write 

<^ltw/o  >  =  |X-aR-aA*A,fJ  dr  dpW(T)W(r  +  p) 

X  <  exp  [  X  (  r  ,  0  )  +  l  (r  +  o,0)]>  .  (2.  78) 

By  using  Eq.  (54)  we  can  rewrite  this  result  as 

(  r  #0  >  =  i  X-3  R-8  A*  A  ^  dT  d^  W  (  r  )  W  (T  +  o  ) 

X  exp  [  -Cx  (0,  0)  +  CJL  (7,  0)  ]  .  (2.79) 

As  noted  in  the  discussion  following  Eq.  (58),  the  log -amplitude  covariance 

function,  (  p,  0  ),  is  in  general  quite  short  range  - much  shorter  range 

than  most  laser  transmitter  diameters  of  interest.  Thus  over  most  of  the 
range  of  the  p  -  integration  the  C^(p,0)  term  can  be  set  equal  to  zero.  This 
leads  to  the  approximate  result  that 
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<  S’ltw/o  >  ~  |X-2  R“a  A*  A  exp  [-C/(0,0)]  J*fd7d£  W(r)  W(r  +  p). 

(2.  80) 

Changing  the  variable  of  integration  from  p  back  to  r  *,  and  then  noting 
that  the  double  integral  then  separates  into  the  product  of  two  identical 
integrals,  each  equal  to  frrDS,  we  get 

<<?ut^o  >c«|x~a  R"3  A*A  exp  [ -Cx  (0,0)]  (^ttD2)2  .  (2.81) 

When  we  substitute  this  result  together  with  Eq.  (74)  into  Eq.  (69), 
and  take  note  of  Eq.  (32),  we  get  the  result  that 

<  Gltm/0  >  «*1tt  (D/X)2  exp  [-C^O.OJ] 

~G0L  exp  [  -Cjj  <0, 0)  ]  ,  (2.82) 

as  the  ensemble  average  adaptive  optics  laser  transmitter  antenna  gain  when 
the  adaptive  optics  does  not  provide  a  random  apodization  correction,  and 
there  is  no  angular  separation  between  the  directions  to  the  target  aim -point 
and  to  the  reference  beacon. 

It  is  obvious  from  consideration  of  Eq.  (82)  that  the  failure  to  provide 
for  the  random  apodization  correction  reduces  the  adaptive  optics  laser  trans¬ 
mitter  antenna  gain  by  a  factor  of  exp  [  -C^(0,  0)].  If  the  log  -amplitude  variance, 
C ^  (0,0),  is  small  enough  this  can  be  ignored,  but  for  nontrivial  values  of  the 
log -amplitude  variance  it  would  be  desirable  to  have  an  adaptive  optics  that 
could  make  random  apodization  corrections - if  possible. 

Before  we  leave  this  subsection  to  take  up  the  case  of  a  laser  trans¬ 
mitter  with  adaptive  optics  capable  of  producing  random  apodization  correction, 
we  wish  to  write  down  the  relevant  result  when  we  do  not  make  the  short  range 
of  the  log -amplitude  covariance  function  approximation  that  we  used  in  going 
from  Eq.  (79)  to  Eq.  (80).  The  r -integration  in  Eq.  (79)  can  be  carried  out 
(inside  the  p  -integration)  to  give  a  factor  of  f  tt  D2  K  (  p  ),  as  noted  in  the  dis  - 
cussion  just  preceeding  Eq.  (21).  Accordingly  Eq.  (79)  can  be  written  as 
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(2.  83) 


(^LTW/O  >  =  R"2  A*  A  (inD2)  f  dp  K  (  p) 

X  exp  [  -Cji(0,  0)  +  C  ^  (p,  0)  ] 

Substituting  this  together  with  Eq.  (74)  into  Eq.  (69)  we  get  the  more  exact 
result  that 

<  GLTIi,0  >  =  X"2  J*  dp  K  (p)  exp  [-^(0,0)  +  C^p,  0)]  .  (2.84) 

Taking  note  of  the  fact  that 

r  dp  K  (  p)  =  j tt D2  ,  (2.85) 

and  using  Eq.  (32),  we  can  rewrite  Eq.  (84)  in  the  very  convenient  form 

<Glt„„  >  exp  r  -C, (0.0)1  .r<£K(<,U«PrCf(?.0>] 

r dp  K  (p) 

=  G0L  exp  (  -C  ,  (0.  0)  ifj  dP  K  (  P)  ^Xp  [C*  j P  ■  °  l  ]| 

L  j'dp  K  (  p)  J.  (2.86) 

To  the  extent  that  the  log -amplitude  covariance  is  very  short  range  com¬ 
pared  to  the  laser  transmitter  aperture  diameter,  D,  the  quantity  in  the  curly 

brackets  will  be  very  nearly  equal  to  unity -  thus  justifying  the  approximate 

result  of  Eq.  (82). 

With  this  result  in  hand  we  are  now  ready  to  take  up  the  case  of  the 
adaptive  optic-  laser  transmitter  with  adaptive  optics  that  can  provide  random 
apodization  correction.  This  is  treated  in  the  next  subsection. 

1.2.  2.4  Ideal  Laser  Transmitter  Adaptive  Optics  With  Random  Apodization 
Co  rrection 

.-TTor  the  case  of  an  adaptive  optics  laser  transmitter  with  the  ability 
to  correct  for  random  apodization  Eq.  's  (65),  (66),  and  (67)  remain  applicable. 
We  rewrite  these  with  the  subscript  LTW  as 
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<  GLTW  >  =  (  <  <?LTH  >  /  <  PliK  >  )  RS 


(2.87) 


<  f?LTH  )  =  2  (  u  L  T  W  ULT  W  )  >  (2.88) 

and 

ul  t  w  =  '  ^  Idr  ULtw  (r)exp{i[kL(r,OT)  +  t(r,0I)]]  ,  (2.89) 

corresponding  to  Eq.  's  (69),  (70),and  (71).  We  also  note  that  corresponding 
to  Eq.  (73)  we  can  write 

(  PLTW  )  =  J  dr  i  <  U*LTW  (7)  ULTM  (r)  >  .  (2.90) 

The  careful  reader  will  have  noted  that  we  seem  to  have  skipped  over  the 
task  of  writing  down  the  expression  corresponding  to  Eq.  (68).  It  would  be 
this  expression  that  would  indicate  the  nature  of  the  random  apodization 
correction.  As  we  shall  see  this  matter  requires  some  careful  consideration 
and  the  appropriate  result  is  not  at  first  apparent. 

Off  hand  it  would  seem  that  the  appropriate  correction  would  be  just 

as  in  the  imaging  case,  as  presented  in  Eq.  (6 1 )  - namely  the  random 

apodization  correction  should  cause  the  the  pupil  function  to  vary  inversely 
in  power  density  to  the  power  density  variation  of  the  beacon  signal  {  In  Eq.  (61) 
the  factor  of  exp  [-jt(r,  0B)  ]  is  the  inverse  of  the  beacon  signal's  power 
density,  which  is  proportional  to  exp  (r,  0B  )  J  } .  If  the  pupil  function,  UL  T  w  (  r  ) 
were  formed  in  this  way  then  the  power  density  variation  would  exactly 
cancel  the  it  (  r,  9T)  induced  variation  in  Eq.  (89)  and  we  would  seem  to  be 
forming  a  diffraction  limited  power  density,  ultw,  at  the  aim -point.  However, 
this  is  not  enough  to  insure  diffraction  limited  laser  transmitter  antenna  gain 
performance  since  these  variations  in  the  pupil  function's  power  density  can 
effect  the  total  laser  power  transmitted,  and  antenna  gain  represents  a  ratio 
of  power  density  at  the  target  aim -point  to  the  total  laser  power  transmitted. 

We  therefore  must  choose  the  nature  of  the  random  apodization  correction  with 
the  average  antenna  gain,  (  GLT  H  ),  in  mind  rather  than  with  the  target  aim - 
point  power  density,  (  &LT  w  ),  in  mind. 
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It  is  not  clear  exactly  how  the  problem  of  defining  the  optimum 
random  apodization  correction  should  be  attempted.  It  appears  to  be  a 
fairly  complex  problem  in  the  calculus  of  variations.  Accordingly  we  shall 
attempt  the  simplier  problem  of  finding  the  value  of  3  in  the  random 
apodization  correction  pupil  function  expression 

ULTW(r)  =  A  W  (  r  )  exp  fik[0(  r  )  +  r  •  0T  ]-0  (  r,  0B)  +  3  X(  r,  0B)],  (2.  91) 

To  do  this  we  need  to  apply  Eq.  (91)  to  obtain  an  expression  for  the  ensemble 
average  laser  power  transmitted  as  a  function  of  3,  and  to  obtain  an  expression 
for  the  ensemble  average  laser  power  density  at  the  target  aim -point,  also 
as  a  function  of  3-  From  these  results  an  expression  for  the  ensemble 
average  antenna  gain  as  a  function  of  3  will  follow  immediately  and  choice 
of  an  optimum  value  for  3  will  be  straight  forward. 

When  we  substitute  Eq.  (91)  twice  into  Eq.  (90)  and  carry  out  the 
appropriate  simplifications  we  get  the  result  that  the  average  laser  power 
transmitted  is 

<  PLth  )  =  iA  A  f  dr  W(r)<  exp  [2  3-M  r,  0B  )  ]  )  .  (2.  92) 

Making  use  of  Eq.  's  (35),  (37),  and  (47)  we  can  write 

<  exp  [2  B  l  (  r,  0B)  ]  >  =  exp  (2  3  l)  exp  [2  3 3  (  (  A  -  i)3  >] 

=  exp  [(2  B8-  23)  <(  l-~if  >] 

=  exp[(2£?  -  23  )C£(0,0)  ]  .  (2.93) 

By  virtue  of  the  stationarity  of  the  statistics  of  the  log -amplitude  this  is 
obviously  not  a  function  of  aperture  position,  r.  Accordingly,  when  Eq.  (93) 
is  substituted  into  Eq.  (92)  the  r -integration  reduces  to  the  evaluation  of  the 
area  of  a  circle  of  diameter  D,  so  that  we  get 

<PLTM>  =  i  A* A  (i ttD2  )  exp  [(  2  £?  -2  3)  CjL  (0,0)]  ,  (2.94) 

for  the  ensemble  average  transmitted  laser  power. 
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To  evaluate  the  laser  power  density  at  the  target  aim -point, 

(  PLT  M  ),  we  first  substitute  Eq.  (91)  into  Eq.  (89).  Appropriately 
simplifying  we  get  the  expression 

uLTW  =  -  ~  exp  [ikR(l+#0T2  )]  A  J’  dT  W  (  r  ) 

X  exp  {iL0  (r,eT)  -0(r,eB)]+f(r,0T)  +  M(r,0e  p.  (2-95) 

At  this  point  we  recall  that  we  are  considering  the  idealized  case  where 
there  is  no  angular  separation  between  the  directions  to  the  target  aim -point 
and  to  the  reference  beacon,  both  being  on  the  z-axis,  so  that  0T  =  0B  =  0. 

This  allows  us  to  reduce  Eq.  (95)  to  the  simplified  form 

u  l  t  h  =  -  exp  ( i  k  R  )  A  fdr  W(r)  exp  [  ( 1  +  3  )  l  (  r,  0  )  ]  .  (2.  96) 

When  we  substitute  this  result  into  Eq.  (88)  and  suitably  simplify,  we  get 
for  the  ensemble  average  laser  power  density  at  the  target  aim -point 

<0LTM  >  =  |  a'"a  X~2  R”2  ff  d?  dr'  W(7)  W(7') 

X  <  exp  f(l  +  3)U(r,  0)  +  l  (r',  0)  }  >  .  (2.97) 

Once  again  making  use  of  Eq.  's  (35),  (37),  and  (47),  we  can  now  write 

<  exp  {  (1 +8  )[  f(r,  0  )  +  f(r',  0)  ]  }  >  =  exp  [2  (  1  +  8  )  f  ] 

x  exp  [  (1  +  B)3[C/(0,0)  +  C£(r'-'?,0)]  } 

=  exp  [  ( 82  -  1  );CX  (0,0)  +  (l +B)S  Cx(r'-r,0)  ]  .  (2.98) 

When  we  substitute  this  result  into  Eq.  (97)  and  replace  the  variable  of 

— >  j  — >  — > 

integration  r  with  p,  [with  p  as  defined  by  Eq.  (72)  J,  we  get 

<  )  =  g  A  A  X'3  R-2  JJ  dr  dp  W  ( r  )  W  (  r  +  p  ) 

X  exp[  (32  -1)CX(0,0)  +  (1  +B)3  Cx(ft0)]  -  (2.99) 

Introducing  the  approximation  that  since  the  log -amplitude  covariance  function, 

C^p.O),  is  a  short  range  function  in  terms  of  the  aperture  diameters  of 

■ — > 

interest  to  us,  we  can  consider  C ^  (  p,  0  )  to  have  zero  value  over  most  of  the 
range  of  integration,  we  can  write 
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<0LTW  )  **%a"a  X-8  R-2  exp[(08-l)  CL  (0,  0)]  l‘J  dr  W  (r  )  W  (7  +7) 

-■=»  i  a" a  r 2  r"s  exp  r ( e2  - i  )Ci  (o,  o)]  rpd7  dr'  w(7)  w  (?') 

«  4  A*A  X-8  R~b  exp  [  (0s  -1)  C£  (0,  0)]  (  £tt  D8)8  . 

(2.  100) 

Obviously,  either  0  =  +  1  or  0  =  -  1  will  result  in  an  ensemble  average 
laser  power  density  of  the  target  aim -point  which  is  independent  of  the 
turbulence  induced  log -amplitude  variance.  But,  the  ensemble  average 
laser  power  transmitted  does  depend  on  0,  as  we  can  see  from  Eq.  (94). 

The  ensemble  average  antenna  gain,  which  is  what  we  want  to  optimize  by 
our  choice  of  0,  is  obtained  by  substituting  Eq.  's  (94)  and  (100)  into  Eq.  (87). 
This  yields  the  result  that 

<  Guth  )  4  tt  (D  /  X)8  exp  [(  -  02  +  2  0  -  1  )C/(0,0)] 

~Gol  exp  [  -  (0  -1  )SCjj  (0,  0)  3  .  (2.101) 

It  is  obvious  from  this  expression  that  the  ensemble  average  antenna  gain, 

(  Glt  h  ),  will  be  maximized  [and  that  this  maximum  will  be  equal  to  the 
diffraction  limited  value  of  antenna  gain,  GD  L  ,  as  defined  by  Eq.  (32)]  when 
0  equals  unity. 

If  we  go  back  to  Eq.  (99)  and  do  not  make  the  approximation  that  the 
log -amplitude  covariance  is  a  short  range  function  relation  to  the  aperture 
diameters  of  interest  then  using  Eq.  (21)  to  allow  the  r -integration  to  be 
performed,  we  get 

(<?LTW>  =  bA'A  X“2  R-8  UttD2)  exp[(fP  -  1  )  Cj  (0,  0)  ] 

x  |  dt  K  (p  )  exp[  (1  +0  f  CiU,  0)  ]  .  (2.102) 

Substituting  this  together  with  Eq.  (94)  into  Eq.  (87)  we  get  as  a  more  exact 
expression  for  the  ensemble  average  antenna  gain 


(  Glt  w  ) 


X"a  exp  [  -(  0  -  1  ) 8  Cx(0,  0)  ] 

x  J  dpK  (  o  )  exp  [(1  +  0)8  (  p,  0  )  ] 
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(2.  103) 


Making  use  of  Eq.  's  (32)  and  (8  5)  this  result  can  be  cast  in  the  convenient  form 


<GLTW  )  =  Gol  exp  [  -(  8-  1  )S  (0,  0)]j 


f‘dpK  (p)  exp  [(1  +  8)2  Ci(p,  0)] 


r  dp  K  (  p) 


(2.  104) 


The  surprising  thing  to  be  noted  from  this  expression  for  the  ensemble 
average  antenna  gain  is  that  while  for  large  enough  aperture  diameters  the 
quantity  in  the  curly  brackets  is  very  nearly  unity  so  that  the  optimum 
choice  of  value  of  8  is  unity  and  the  maximum  possible  value  of  antenna 
gain  is  the  diffraction  limited  value,  for  small  enough  values  of  the  aperture 
diameter  the  situation  is  quite  different.  For  very  small  values  of  aperture 
diameter  the  ensemble  average  antenna  gain  becomes  very  nearly  equal  to 
Gol  exp  [  4  8  Cjj  (0,  0)  ].  In  this  case  very  large  values  of  8  appear  desirable 
and  the  ensemble  average  antenna  gain  can  be  greater  than  the  diffraction 
limited  value.  However,  this  is  a  matter  of  having  the  laser  transmitter 
emit  the  largest  amount  of  power  at  those  instances  when  the  power  density 
of  the  beacon  signal  indicates  that  the  "connectivity"  between  the  laser  trans¬ 
mitter  aperture  and  the  target  aim -point  provided  by  the  turbulence  is  greatest. 
Though  it  seems  interesting  to  note  here  this  somewhat  surprising  possibility, 
we  believe  that  for  a  number  of  practical  reasons  the  possibility  of  exploiting 
this  feature  is  relatively  uninteresting. 


Before  leaving  the  subject  area  of  laser  transmitter  adaptive  optics 
with  random  apodization  correction  capability  it  is  perhaps  worth  remarking 
on  the  implications  of  having  p  =  1  in  Eq.  (91).  This  means  that  we  have 

ULTW(r)  =  A  W  (  r  )  exp  fi  k  [Q(  r  )  +  r  •  ST]  -i0  (  r,  06)  +  L  (  r  0B)  1.  (2.105) 


Making  use  of  Eq.  (34)  we  can  rewrite  this  as 

ULTW(r)  =  A  W  (  r  )  exp  fik  [o(  r  )  +  r  •  0T  ]  -i  *  (  r,  0B  )  1  .  (2.106) 

Comparing  this  with  Eq.  (62)  for  the  imaging  system  with  random  apodization 
compensation,  we  see  that  while  the  adaptive  optics  of  the  imaging  system 
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results  in  subtraction  of  the  complex  phase  of  the  beacon,  for  the  laser 
transmitter  it  is  the  conjugate  of  the  complex  phase  of  the  beacon  that 
is  subtracted. 


With  all  of  these  results  in  hand  for  imaging  system  and  laser 
transmitter  performance  under  the  ideal  conditions  of  no  angular  separation 
between  the  directions  from  the  aperture  to  the  target  and  to  the  beacon, 
with  both  target  and  beacon  lying  on  the  z  -axis,  we  are  now  ready  to  turn 
our  attention  to  the  cases  of  basic  interest,  when  the  angular  separation 
is  non -zero  and  there  is  a  possibility  of  anisoplanatism  effects.  However, 
before  we  can  evaluate  system  performance  in  such  a  situation,  we  must 
first  develop  expressions  for  the  joint  statistics  for  propagation  in  two 
distinct  directions.  This  is  taken  up  in  the  next  section. 


I 


The  analytic  foundation  for  the  propagation  theory  that  we  shall 
develop  in  this  section  will  be  found  in  an  earlier  paper  of  oursb  titled 
"Spectral  and  Angular  Covariance  of  Scintillation  for  Propagation  in  a 
Randomly  Inhomogeneous  Medium"-  In  the  following  we  shall  refer  to  this 
paper  as  SAC.  The  portion  of  SAC  concerned  with  the  angular  covariance 
is  directly  relevant  to  our  requirements  here,  though  it  needs  some  exten  - 
sions  and  modifications  to  provide  the  propagation  theory  basis  for  the 
work  undertaken  here.  (There  is,  in  addition  a  requirement  for  one  non¬ 
trivial  correction.  )  While  the  results  of  SAC  are  given  for  the  log -amplitude 
and  for  the  phase  covariance,  no  results  are  given  for  the  cross -statistic s 
between  phase  and  log -amplitude.  We  shall  need  this  latter  result  to  estab¬ 
lish  the  isotropy  of  the  cross -statistics - a  fact  which  we  shall  use  later 

and  which  is  not  apparent  from  any  generalized  reasoning.  In  addition,  while 
the  results  of  SAC  are  developed  for  infinite  plane  wave  propagation,  to  allow 
us  to  apply  our  work  to  the  case  in  which  the  adaptive  optics  system  and  the 
target/beacon  are  within  the  atmosphere,  we  shall  need  corresponding  results 
for  spherical  wave  propagation  theory.  Our  approach  to  developing  the  spheri¬ 
cal  wave  results  will  be  heuristic  arguments  extending  the  infinite  plane  wave 
results.  It  will  be  the  detailed  development  of  the  infinite  plane  wave  results 
that  we  shall  concentrate  on,  and  only  at  the  very  end  of  this  section  will  we 
introduce  the  heuristically  justified  modifications  that  will  make  these  results 
applicable  to  the  spherical  wave  case. 

In  what  follows  we  shall  not  repeat  the  portion  of  SAC  concerned  with 
the  equations  for  propagation  through  a  random  medium  and  the  development 
of  the  random  wavefunction.  We  shall  simply  quote  these  results  and  then 
apply  them  to  the  development  of  the  statistics.  We  shall,  however, go  over 
most  of  the  notation  as  we  wish  to  make  some  (simplifying)  changes  in  the 
notation.  We  start,  in  the  next  subsection  with  a  presentation  of  the  notation. 


1.3.1  Notation  and  Random  Propagation  Results  from  SAC 


The  basic  random  function  of  interest  to  us  is  the  previously 
introduced  complex  phase,  t  (  r,  0  ),  with  its  real  and  imaginary  parts, 

0  (7,  7 ),  the  real  phase,  and  L  (r,7),  the  log -amplitude  variance,  which 
by  copying  Eq.  (2.  34)  we  write  here  as 

ijt  (7,  7)  =  0  (  r,  0  )  -  i  JL  (  r,  0  )  .  (3.  1) 

Here  r  denotes  a  position  on  the  aperture  plane  at  which  the  complex  phase 

is  measured,  and  0  indicates  the  direction  from  the  center  of  the  aperture 

to  a  monochromatic  point  source  of  wavelength  X  at  a  z-axis  distance  equal 

to  R.  The  angle  0  is  a  measure  of  the  deviation  of  the  direction  from  the 

z-axis,  which  axis  is  perpendicular  to  the  aperture  plane.  (For  a  point 

source  on  the  z-axis  the  angle  0  is  equal  to  zero.  )  In  addition  to  the  complex 

phase,  ty,  the  other  basic  random  function  which  we  have  to  consider  is 

n  ( 7 ,  z),  the  random  deviation  of  the  refractive -index  from  its  nominal  mean 

value  of  unity.  (It  is  this  random  function  that  we  are  referring  to  when  we 

speak  of  turbulence.  )  The  refractive -index  denoted  by  n  (  r,  z  )  is  measured 

at  a  distance  z  from  the  aperture  plane  along  the  z  -axis,  and  at  a  position 
— — ► 

r  displaced  from  the  z-axis. 

In  Eq.  (2.  47)  we  introduced  the  covariance  of  the  log  -amplitude,  . 

We  repeat  that  definition  here  in  a  slightly  modified  form  along  with  definitions 
of  all  the  other  statistical  function  we  shall  be  concerned  with.  These  are  as 
follows;  the  log -amplitude  mean  value 

7  =  <  i  (7,7)  >  ,  (3. 2) 

the  real  phase  mean  value 

J  =  <  0  (7,  7  )  >  ,  (3.  3) 


the  log -amplitude  convariance 


the  phase  covariance 


C#  (  P»  «?  )  =  ([0{r  +  p,  §«?)-0][  0(r,-|i?)-0]>,  (3.  5) 

the  log -amplitude:  phase  cross -covariances 

CX(p,?)  =  <  [jMr+7,  0) -l][0(r,  -id)  -0  j>  ,  (3.6) 

and  the  ref ractive -index  covariance 

CN  (p,  §)  =  (  n  (  r  +  p,  z  +  5  )  n  (  r,  z  )  >  .  (3.7) 

We  shall  have  occasion  to  work  with  the  Fourier  transforms  of  the  above 
four  covariance  functions.  For  the  log -amplitude  covariance,  phase  co- 
variance,  and  log -amplitude:phase  cross -covariance  these  are 

(  rr,  i?  )  =  J*  dp  (p.tf)  exp  (-i  rr  •  p  )  ,  (3.  8) 

F#  (a,*)  =  j  dp  C#  (p,  j?)  exp  ( -  i  a  ’  p  )  ,  (3-9) 

and 

Fjj*  (  a,  )  =  J'  dp  (  p,  0)  exp  ( - i  a  *  p  >  ,  (3.  10) 

respectively.  For  the  refractive -index  covariance  we  shall  consider  both 
a  two-dimensional  Fourier  transform,  namely 

FN  ( <r,  ?  )  =  J  dp  CN  (  p,  ?  )  exp  (  -  i  o  •  p  )  ,  (3.  11) 

and  a  three-dimensional  transform 

3n  (  c,  s  )  =  fj  dp  d§  CN  (  p,  §  )  exp  l  -  i  (  <?  '  P  1  s  §  )  ]  .  (3.  12) 

It  is  immediately  obvious  from  consideration  of  Eq.  's  (11)  and  (12)  that 

3n  (  '7,  s  )  =  /  d§  FN  (jt,  f  )  exp  (  -i  s  5  )  .  (3.  13) 

As  noted  in  SAC  experimental7  and  theoretical  evidence6*  9  based  on  the 
K olmogoroff 1  s  theory  of  turbulence  in  the  inertial  subrange  and  Tatarski's 
arguments  of  a  conserved  passive  additive  nature  for  the  optical  aspects  of 
turbulence  lead  to  the  conclusion  that 

On  (  ">  s  )  =  8.  16  CNS  (rr8  +  s3  f 1 1/6  ,  (3.14) 


where  CN2  is  a  scalar  quantity  called  the  refractive -index  structure  constant 
which  serves  as  a  measure  of  the  local  optical  strength  of  turbulence.  (The 
dimensionality  of  Cm2  is  length  to  the  minus  two -thirds  power.  ) 

In  addition  to  these  Fourier  transforms  of  statistical  functions  SAC 
introduces,  we  shall  have  to  make  use  of  Fourier  transforms  of  certain 
random  functions.  These  Fourier  transforms  are  as  follows:  for  the 
complex  phase  we  have 

'Met,  0)  =  f  dr  i|i(r,«)  exp(-ir*r)  ,  (3.15) 

for  the  log -amplitude  we  have 

L(rr,  9)  =  f  dr  X  (  r,  0  )  exp  (  -  i  ^  •  r  )  ,  (3.16) 

for  the  real  phase  we  have 

$  (  cr,  3  )  =  |dr  0  (r,9)  exp(-irr*r)  ,  (3.17) 

and  for  the  refractive -index  we  have 

N  (  cr,  z  )  =  J  dr  n(r,  z)  exp  (  -i  c  •  r  ) 

It  follows  directly  from  Eq.  's  (1),  (15),  that 

L(a,?)  =  £i[  7(o,0)  -  7*( -a,  9)  ] 
and  from  Eq.  's  (1),  (15),  and  (17),  that 

Ma,0)  =  i[’f(o.9)  +  ’f*(-o,0)]  .  (3.20) 

The  solution  of  the  propagation  equation  carried  out  in  SAC  leads 
* 

to  the  result  that 

*  Eq.  (21)  above  is  in  some  significant  ways  different  from  Eq.  (20)  in  SAC 
from  which  it  is  taken.  Some  of  these  differences  are  due  to  the  fact  that 
7  as  given  by  Eq.  (1)  here  is  equal  to  7  as  defined  in  SAC  times  i.  A 
second  cause  of  difference  is  that  while  in  SAC  the  integration  runs  from 
the  source  to  the  plane  where  7  is  measured  (with  the  variable  of  inte¬ 
gration  being  equal  to  zero  at  the  source  and  equal  to  R  at  the  measure¬ 
ment  plane),  here  the  integration  runs  from  zero  at  the  (aperture)  plane 
where  7  is  measured  to  z  =  R  at  the  source. 


(3.  18) 

(3.  19) 
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where 


8 

2  ik  j  dz  N  (  rr,  z)  exp  (  -ikz  )  sin  [y  (  a,  8  )  z  ],  (3.  21) 
o 

k  -  2  TT  -  a  •  9  -  I  k  02  .  (3.  22) 

Before  leaving  this  subsection  it  will  be  convenient  to  rewrite  Eq.  (21)  as 

Y  (  o,  0  )  =  -k  j*  dz  N  (rr,  z  )fexp  {i  [y  (  o,  9  )  -k]z  }-exp  {  -  i[v  (  o,  9)  +  k]z]  .. 

o'-  / 

(3.  23) 

and  from  this  to  write 

(-rt,  9  )  =  -k  r  dz  N*  (-rr,  z)  (exp  {  -i[  Y(-rr,  0)-k]  z  }  -exp  {i  [v  ( -u,  8)  +k  ]  z  ]\ 
o  '  ' 

(3.  24) 

By  virtue  of  the  fact  that  the  refractive -index  variation  is  real  (so  that 
n  =  n  )  it  follows  from  Eq.  (18)  that 

N*  (-",  z)  =  N  (~,  z)  ,  (3.  25) 

so  that  we  can  rewrite  Eq.  (24)  as 

Y  ( 9  )  =  -k  ['  dz  N  (  <7,  z  )  (expf  -i  [  v  (  -o,  8)  -k  ]z  }  -exp  {i  [v(  -o»  6)  +  k  ]z  1  , 

(3.  2  6) 

When  Eq.  ‘s  (23)  and  (26)  are  substituted  into  Eq.  's  (19)  and  (20)  and  the 
obvious  simplifications  are  made,  we  get  the  results  that 

p 

L  (  9  )  =  -  |  i  k  J  dz  N  (  t,  z  )  ^  exp  { i  [  y  (  rr,  9  )  -k]z  1 

o 

-  exp  {  -i  [  V  {  r-,  9  )  +  k  1  z  1  -  exp  f-i[v(-o,9)-k]z} 

+  exp  {  i  [Y  (  -  rr,  8  )  +  k  ]  z  1  j  ,  (3.27) 

and 
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Hf  (n,9)  = 

V  (0,0)  = 


1 


_  _  »  _  _ 

$  (  a,  0  )  =  -|ik  J*  dzN(«,z)  ^exp  [i  [v  (  rr,  0  )  -k  ]  z  } 

o 

-  exp  {  -  i  [  ^  ( <r,  0)+k3z}  +  exp  {  -i  [  V  (  -n,  0  )  -  k  ]  z  ] 

-  exp  f  i  [Y  (-CJ,  0  )  +k  ]  z  }  ^  .  (3.28) 

At  this  point  we  note  that  of  the  four  exponentials  in  each  of  the 

above  two  integrands,  as  we  can  see  from  consideration  of  Eq.  (22),  and 

— ♦ 

of  the  facts  that  k  is  much  greater  than  cr  and  that  the  magnitude  of  0 
is  very  small,  two  of  the  exponentials  oscillate  very  rapidly  (at  a  "rate"  of 
about  2  k),  while  the  other  two  oscillate  much  more  slowly.  Such  rapid 
oscillations  can  not  contribute  to  the  value  of  the  integrals  and  accordingly 
we  can  make  the  approximations  that 

L(a,^)=  -  5  ik  fdzN(fr,  z)  ^exp  {i  [v  (  o,  0,  )-k]z} 

0 

-  exp  {  -i  [V  ( -  n,  0  )  -k  ]  z  }  j  ,  (3.  29) 

and 

p 

$  (  a,  0  )  a  -|k  f  dz  N(  o,  z)  ^exp  {i  [V  (  cr,  0  )  -k  ]  z  } 

0 

+  exp  {  -i  ['f  ( -n,  0  )  -k  ]  z  }  y  .  (3.  30) 

Making  use  of  Eq.  (22)  we  can  recast  these  results  in  the  form 

P 

L  ( er,  0  )  =  -  4  i  k  f  dz  N  (  a,  z  )  ( exp  {  [  -i  (|  k  02  +  £  ns  /k)  -  i  rr  •  0  ]  z  ] 
o 

-  exp  {  [i(|k0s  +|rrs/k)  -i^- 0  ]  z  }  /  ,  (3.31) 


p 

$  ( cr,  0  )  =  -  |k  T  dzN(c.z)  (exp  {[-i(£k02+i  a2/k)-i  cr  •  0  ]  z  1 
0 

+  exp  {  [i  (|k  02  +  i  rj2  /k)  -i  c  •  9  ]  z  }  j  .  (3.  32) 

This  is  as  far  as  we  can  go  with  the  analysis  of  the  random  propagation 
functions.  To  proceed  beyond  this  point  we  have  to  start  working  with  the 
statistical  functions.  To  do  this  we  need  to  make  use. of  a  special  theorem 
developed  in  SAC.  This  is  presented  in  the  next  subsection. 
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1.3.2  Fourier  Transform  Statistics 

In  this  brief  subsection  we  wish  to  prove  a  theorem  about  the 
ensemble  average  of  the  product  of  Fourier  transforms  of  random 
functions.  We  will  show  that  an  integral  over  such  an  average  gives 
rise  to  a  power  spectrum  (or  cospectrum)  for  the  (two)  random  functions, 
with  a  Dirac  delta  function  like  property  associated  with  the  average. 

• 

We  consider  two  real,  stationary  random  functions,  f  (x;o)  and 
g  (x;3),  where  x  is  an N -dimensional  variable  and  a  and  3  are  some 
multicomponent  parameter  "vectors".  (The  functions  f  and  g  need  not 
be  distinct,  but  for  generality  we  allow  them  to  be  so.  )  We  define  the 
cros s -covariance  as 

CFG  (x'-x  ;  a,  3)  =  <  [f  (x';  a) -f  (ct)  ]  [g  (x  ;  3) -g  (3)  ]  >  ,  (3.33) 

where 

f"(a)  =  <f(x;a)>  ,  (3.34) 

g<3)  =  <  g  ( x ;  3 )  >  .  (3.35) 

In  addition  we  write  the  Fourier  transforms  of  the  two  random  functions  as 

F  (y,  a)  =  P  dx  [f  (x;  a)  -f  (a)  ]  exp  (-ix  •  y  )  ,  (3.36) 

G(y,  3)  =  [  dx  [g  (x;  3)  -g  (3)  ]  exp  (  -ix  *y  )  .  (3.37) 

We  also  note  that  we  can  write  the  cospectrum  associated  with  these  two 
random  functions  as 

rFG  (y  ;  a,  3)  =  P  dx  CFG  (x  ;  a,  3)  exp  (  -ix  •  y  )  .  (3.38) 

We  now  introduce  the  function  h  (y)  to  represent  any  nonrandom,  reasonably 
well  behaved  function  and  consider  the  problem  of  evaluating  the  expression 

fdy'(F  (y' ;  a)  G  (y  ;  3)  >  h(y') 

To  carry  out  the  evaluation  we  start  by  making  use  of  Eq.  's  (36) 
and  (37).  This  allows  us  to  write 
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a.  • 


At  this  point  we  introduce  sum  and  difference  coordinates 


p  =  x'  - x  , 

q  =  4  (x'  +  x) 

and  make  these  the  variables  of  integration  instead  of  x  and  x  ,  in  Eq.  (39)* 
Doing  this  and  making  use  of  Eq.  (33)  we  can  now  write 

f  dy'  <  F*(y';a)  G(y;|)  >  h  (y')  =  fff  dp  dq  dy' 

xC,G  (p  ;  a,  P)h(y')  exp  [ip  •  )  +  iq  '  (y'*y)  1  .  (3.  40) 

where  the  exponential's  argument  has  been  rearranged  according  to  the 
formula 

ab  -  cd  =  (a  +  c)  )  +  ^  +  *  (3-41) 

Now  we  take  note  of  the  fact  that  according  to  Eq.  (2.  11)  we  can  write 

J*J  dq  dy'[h(y')  exp  (ip -y 72)]  exp[iq  •  (y'  -y)] 

=  (2  tt)n  h  (y)  exp  (ip  •  y/2)  ,  (3.42) 

where  N,  we  recall,  is  the  dimensionality  of  our  variables.  When  we 
substitute  Eq.  (42)  into  Eq.  (40)  and  simplify  as  appropriate  we  obtain 
the  result  that 

/  dy'  (  F*(  y';  a  )  G  (  y  ;  6  )  >  h  ( y ")  =  (2  n  )N  h  (  y  ) 

x  J  dp  C,  6  (p  ;  n  ,  0  )  exp  ( i  p  •  y  )  .  (3.43) 

Making  reference  to  Eq.  (38)  and  noting  that  since  f  and  g  are  real  functions 
then  so  is  Cr  Q  ,  we  see  that  the  integral  on  the  right -hand -side  of  Eq.  (43) 


can  be  identified  with  Tr  G  .  Thus  we  can  write 

fdy'(F  (y';a)G(y;0)  >  h(y')  =  (  2  tt  )n  h  (  y  )  Vr  G  *(  y  ;  a  ,  ? ).  (3.  44) 

This  is  the  theorem  we  wished  to  prove  in  this  subsection. 

Three  alternate  forms  of  this  result  are  as  follows: 

J*d7'<F(7';a)  G*(y;g)>  h(y')  =  (  2  "  )N  h(y)  rw  (yin.?)  ,  (3.45) 
J*dy'<F*(y:a)G(7';e)>  h(y')  =  (2tt)n  h(y)  rF  G  (y;a,  B),  (3.46) 


Pdy  '  <  F  (y  ,  o  )  G  (y  3  )  >  h  ( y  ')  =  (  2  tt  )n  h  (y  )  rFG  (y  ;  a  ,  ?)  .  (3.  47) 

*  -*  -»  -» 

Eq.  (46)  is  obtained  by  noting  that  replacing  F  (y  ;a)G(y;3)by 

— *  .  — .  j  — , 

F  (  y  ;  a  )  G  ( y  ;  (3  )  would  have  altered  the  right  hand  side  of  Eq.  (40)  only 
to  the  extent  of  interchanging  y  and  y  in  the  exponential,  which  as  we  can 
see  from  consideration  of  Eq.  (2.  11)  would  result  in  no  change  at  all  in  Eq.  (44). 
Eq.  (45)  is  obtained  as  the  complex  conjugate  of  Eq.  (44),  while  Eq.  (47)  is 
obtained  in  the  same  way  from  Eq.  (46). 

With  these  results  in  hand  we  are  now  ready  to  develop  statistical 
propagation  results  from  our  random  propagation  expressions.  We  take  this 
up  in  the  next  subsection. 


1.3.3  Statistical  Propagation  Results 

We  can  write  the  log -amplitude  covariance  function  as  the  inverse 
Fourier  transform  of  F^  as  defined  by  Eq.  (8).  We  can  see  that  this  is  to 
be  written  as 


Ci  (  o,  s?  )  =  (  2  tt  )-2  f  d  u  (  a,  d  )  exp  ( i  a  •  p  ) 
Similarly  we  can  write  for  the  phase  covariance 

C«(p,*?)  =  (2tt)"2  f  der  F#  (<?,  £)  exp  (irr  •  p  ) 

and  for  log -amplitude:  phase  cros  s -covariance 

(  P>&)  =  (  2  tt  )“  2  f  d  <t  &  )  exp  (  i  r  •  "p  ) 


(3.  48) 


(3.  49 ( 


(3.  50) 


r 
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Obviously  what  we  need  to  be  able  to  evaluate  these  three  (cross  -covariance 
functions  is  expressions  for  the  corresponding  Fourier  transforms,  and 
these  we  can  develop  by  a  judicious  combination  of  the  results  of  the  two 
preceding  subsections. 

Making  use  of  Eq.  's  (45)  and  (27)  we  can  write 

(2n)3  Ft  (o  ,?)  =  fdo'  <L(?,|J)  L*  (p,  -£*)  >  .  (3.51) 

Now  making  use  of  Eq.  (31)  we  can  recast  this  in  the  form 

F^  (~,  ?)  =  (4k/n)£  If  dz'dz  j*dr'-<  N(p',  z')  N  (  p,  z  )  > 

00 

X  (  exp  {[-L(i-ki9a  +  ^rt^2/k)  -  j  i<t '•  t)  ]  z'] 

-  exp  {[i(ik  £2  +  |r'2/k)  -f  ip'-?]z'}) 

X  (  exp  f[if|-kjJa+|'CT2/k)  - 1  i  p  •  i?  ]  z  1 

-exp  f[-i(i.k02  +tcr2/k)  -tip  -?]z}).  (3.52) 
Proceeding  in  exactly  the  same  way,  we  also  write 

(  2  tt  )2  F,  (p,  ?)  =  J  dp'  <  Ho',  I?)  **  (a,  -fc?)  >,  (3.  53) 

and 

(2n)2  F*#£,  ?)  =  fdp'  <L(p',fc?)  *(p,  -*?)>.  (3.54) 

Now  making  use  of  Eq.  (32)  as  well  as  Eq.  (31)  we  can  recast  these  as 

F*  (  p,  &  )  -  (ik/n)3f'rdz#dzrdcr'(N(rr',  z  ')  N  (  ~,  z  )  > 

00 

X  (  exp  {  [  -i  fg-kt?2  +  £p'2/k  )  -fia''!>]z'] 

+  exp  f  [  i  (%-  k  t?s  +  |  a'2  /k )  -  |ip'-«?]z''}) 

X(  exp  f  [  i  (i-k  t?3  +  ^  n3  /  k)  -  ^  ip  •  ]z  ] 

+  exp  f  [ -i  (j  kt?3  +  ^  p3 /k)  -  ^  i  rr  ■  ]z  1  ),  (3.55) 
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and 


F^#  (  o,  #)  =  (ik/n)2  Jfdz'dz  Jd  o'  <  N  (  n',  z')  N  (  ^,  z  )  ) 

o  o 

X  (  exp  {  [-i(i-k«?a  +|  o'2/k)  -  \  i  a  '  ' «?  ]  z  '  ) 

-  expf  [i(i-k^2  +|a'2/k)  -iia'-?]z'}) 

X  (  exp  { [  L  (-k  k  «?2  +  J  a2/ k)  -  f  in  •  s>]z] 

+  exp  f[-i(±  kd3  +|  io  -?]z  1)  .  (3.56) 

From  consideration  of  Eq.  's  (11)  and  (45)  it  is  easy  to  see  that 

Jda'<N(^z')N*Uz)>E(^)  =  (2n)a  FN  (  a,  z  '  -z  )  E  (o  ),  (3.  57) 

where  E  (  n)  is  any  well  behaved  function.  Making  use  of  this  result  we 
can,  with  trivial  effort  carry  out  the  n  -  integrations  in  Eq.  's  (52),  (55), 
and  (56).  Thus  we  can  write 

F^  (  n,  J?  )  =  (I  k  f  JT  dz  '  dz  Fn  (  o,  z'  -z  )  exp  [  -  |i  a  •  0  (z'  +  z  )  ] 
00 

X  {exp  [- i u(z ' -z) 3  -  exp  [  - i u(z  *  +  z )  ] 

-  exp  [i  u(z  '+  z)  ]  +  exp  [i  u(  z  *  -  z  )  ]  }  ,  (3.58) 


R  R 

F*(o,t?)  =  (|kf  JJ  dz'  dz  FN(  a,  z'-z  )  exp  [-|  i  rr- (z"+z  )  ] 

00 

X  [exp  [-i  u(z  '  -z)  3  +  exp  [  -ip(z'+  z)  ] 

+  exp  [  i  u(z  '+  z)  ]  +  exp  [  i  u(z  '  -  z)  ]  }  ,  (3.  59) 

and 
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F (  a  ,  &  )  =  (r  k)2  I  f  dz  *  dz  F  N  (  re ,  z  *  -z  )  exp  [  -  J  i  o  ■  d  (z  '+  z)  ] 
oo 

x  {  exp  [-i  p  (z  '-z)  ]  +  exp  [-i  p(z  '+  z )  ] 

-  exp  [ip(z '+  z)  ]  -  exp  [  i u(z  *  -z)  ]  } ,  (3.  60) 

where 

p  =  4  k«?a  +  §r3/k  .  (3.  61) 

At  this  point  we  shall  introduce  new  variables  of  integration,  namely 

u  =  z  *  -  z  ,  (3.  62) 

and 

v  =  i  ( z'  +  z  )  .  (3.  63) 

Because  of  the  finite  range  of  the  refractive -index  correlation  function 
(which  range  we  assume  is  much  less  than  the  propagation  path  length  R) 
we  can  consider  the  limits  of  the  u -integration  to  run  from  to  +  oo,  while 
the  limits  of  the  v  integration  run  from  0  to  R-  Thus  we  can  rewrite 
Eq.  's  (58),  (59),  and  (60)  as 

Fjl  (  a,  ;?)  =  (®  k)3  Jdv  Jdu  FN  (  a,  u  )  exp  (  -  i  a  •  v  ) 
o 

X  [  exp  (  -i  pu  )  -  exp  (  -2  i  p  v  )  -  exp  (2  i  p  v  )  +  exp  ( i  uu  )  ]  , 

(3.  64) 

R 

F#  (  £ )  =  (|  k)3  J  dv  Jdu  F»-  (  o,  u  )  exp  (  -i  a  •  0  v  ) 

6 

x  [  exp  (  -i  pu)  +  exp  (  -2  i  pv)  +  exp  (2  i  pv)  +  exp  (i  pu  )  ]  ,  (3.6  5) 

and 

F £ #  ( a,  j?  )  =  (| k)3  J  dv  Jdu  FN  (  c,  u  )  exp  ( -i  a  •  v  ) 
o 

x  [  exp  (  -i  au  )  +  exp  ( -2  i  pv  )  -  exp  (2  i  pv)  -  exp  (ipu)]  .(3.66) 
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Rather  than  carry  the  F ^  -  term  any  further  we  shall  now  show 
that  it  is  identically  equal  to  zero.  This  follows  from  the  fact  that  since 
i  and  t  are  each  real  valued  (random)  quantities,  then  C ^  must  also 
be  real.  From  this  fact,  together  with  Eq.  (10),  we  can  see  that 

Fj it  (  -  c.  *)  -  F/#  (a,  ?)  •  (3.  67) 

Similarly  since  the  random  ref  ractive -index,  n,  and  its  covariance 
function,  CN,  are  real  valued,  it  follows  from  consideration  of  Eq.  (11) 
that 

F?  (  -  a,  l)  =  Fn  (Z,  ?)  .  (3.  68) 

We  note  that  according  to  Eq.  (61)  the  value  of  g  is  unchanged  when  9  is 
— > 

replaced  by  -  n.  This  means  that  we  can  write  as  an  alternative  form  of 
Eq.  (66)  the  equation 

Fjj0  (  -  a  ,  *?  )  =  (|  k)2  f  dv  J  du  Fn'  (  -  n  ,  u  )  exp  (  -  i  a  •  0  v  ) 

o 

X  [  exp  (imu)  +  exp  (2  i  g  v)  -  exp  (  -2  i  u  v)  -  exp  ( -  i  p,u)  ]  . 

(3.  69) 

Making  use  of  Eq.  (68)  we  can  recast  this  expression  in  the  form 

(-?,$)  =  -  (g  k)a  J  dv  f  du  F«(",u)  exp  (  -i  a  •  <9  v  ) 
o 

X  [  exp  (  -i  u  u  )  +  exp  (  -2  i  gv  )  -  exp  (2  i  gv  )  -  exp  (i  uu)  ]  . 

(3.  70) 

It  is  obvious  that  Eq.  's  (66)  and  (70)  are  compatible  with  Eq.  (67)  only  if 

FU  (c,  *)  *  0  (3.  71) 

According  to  Eq.  (50)  this  means  that 

(■?.?)=  0  .  (3.72) 

We  are  now  ready  to  return  to  the  evaluation  of  Fjj  and  F*  and  then  of 
C i  and  C*  . 
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We  can  carry  out  the  u -integration  in  Eq.  's  (64)  and  (65)  by 
noting  that  in  accordance  with  Eq.  (13)  the  u -integrations  can  be  con¬ 
sidered  to  correspond  to  taking  a  Fourier  transformation  along  the  z-axis. 
This  will  convert  the  two  dimensional  Fourier  transform  FN  ( -  ,  u  )  t 
the  three-dimensional  Fourier  transform  (^,  s)  ,  where  s  will  take 
a  value  of  either  0  or  ±2 u-  We  argue  that  u  is  (except  for  some  very 

extreme  conditions  of  no  interest  to  us)  much  smaller  than  the  magnitude 
► 

of  o  .  so  that  in  accordance  with  Eq.  (14)  we  can  consider  the  three -dimen  - 
sional  Fourier  transform  in  all  cases  to  have  arguments  of  ('-,0).  It  is 
convenient  to  write  this  as  ( a)  where 

(o)  s  SH  (7,0)  ,  (3.  73) 

and  according  to  Eq.  (14) 

%  (o)  =  8.  16  C„2  a~1V3  •  (3.  74) 

When  we  carry  out  the  u -integrations  in  Eq.  *s  (64)  and  (65)  we  get 

p 

F^(  o,»)  =  ik2  f  dv  ?N  (a)  exp(-i^.?v)  [1  -cos  (2uv)]  ,  (3.75) 

o 

and 

F*  (  cr ,  t?  )  =  i  k2  f  dv  (  rr)  exp  (-io  • «?  v  )  [1+  cos  (2U  v  )  ]  .  (3.  76) 

o 

Now  substituting  these  results  into  Eq.  1  s  (48)  and  (49)  we  obtain 
the  equation 

Cx  (1>. «?)  =  '^'g1  ~  (k/n)2  J  dv  CN2  f'dr  ?-lV3  [1  -  cos  (2  uv  )  ] 

0 

X  exp  [  i  o-(  p  -  i?  v  )  ]  ,  (3.77) 

and 
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c# (p, *)  = 


8.  16 


8  (k/rr)2  r  dv  CN3  J*da  o  -1 3/3  [1  +  cos  (2  )  ] 

0 

(3.  78) 


X  exp  [in-(p-t?v)3 
Making  use  of  the  well  known  formula  that 

CD 

r  d  a  g  (  rr)  exp  (  i  a  ~  •  x  )  =  2n  r  a  d  CT  g  (  o)  J0  (  a  rr  x  )  ,  (3.  79) 


where  g  (  n)  is  any  function,  then  we  can  rewrite  Eq.  ‘s  (77)  and  (78)  as 

,  R  CD 

Ci  (p,$)  =  —2.  k2  Jdv  cN2  fda  n~e'3  Jo  (a  I  P  -*  v|) 


x  [1  -cos(2uv)]  ,  (3.80) 

and 

C,  (7.^)  =  k2  "dv  C/  Td-  rr'8/3  J0  (-  |  o’-?v|  ) 

X  [1  +  cos  (2p,v)]  .  (3.  81) 

These  two  equations,  together  with  Eq.  (72),  represent  our  basic  results 
for  infinite  plane  wave  propagation.  In  the  next  subsection  we  shall  present 
heuristic  arguments  for  the  changes  needed  to  make  these  results  apply  for 
point  source/spherical  wave  propagation. 


1.3.4  Point  Source/Spherical  Wave  Propagation 

Without  going  through  an  even  more  extensive  derivation  than  that 
of  the  preceding  subsection  we  can  infer  the  appropriate  spherical  wave 
results  from  our  infinite  plane  wave  results.  To  do  this  we  start  with  Eq.  's 
(80)  and  (81)  [as  well  as  Eq.  (72)]  for  the  infinite  plane  wave  case,  and 
apply  the  following  heuristic  argument.  The  significant  difference  between 

the  infinite  plane  wave  case  and  the  spherical  wave  case  is  that  for  the 

- — ► 

infinite  plane  wave  case  o  not  only  denotes  the  separation  of  two  points  on 
the  measurement  (aperture)  plane  but  also  denotes  the  separation  of  two  rays 
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traveling  to  those  two  points  quite  independent  of  where  along  the 

v -integration  (propagation)  path  we  examine  the  matter.  On  the  other 

hand  for  the  spherical  wave  case  the  separation  between  the  two  rays 

is  p  (R-v)/R,  which  varies  with  v  along  the  integration  (propagation) 

■  ■  ► 

path.  We  argue  that  the  presence  of  p  inside  the  v -integrand  in  our 
infinite  plane  wave  results  is  supposed  to  correspond  to  the  separation 
of  the  two  rays  leading  to  the  measurement  points  as  a  function  of  v. 
Accordingly  we  would  write  for  the  spherical  wave  case 

CX(P.3>  =  IV^Jdv  C,2  Jdr  o'*3  Jo  (a  17  ~  -  &  I  } 

o  o 

X  [1  -cos  (2uv)  ]  ,  (3.  82) 

and 


C,  (pj)  =  ~ ^  k2  fdv  Q,2  |dff  a'8/3  J0  (a  |p  ~  |  \ 

o  o 

x[l+cos(2p,v)]  .  (3.83) 

And  of  course,  we  still  have  the  result  that 

C U  (p,  £)  =  0  .  (3.  84) 

It  is  convenient  to  introduce  the  notation 

F^v  =  ~p  (R-v)/R  ,  (3.  85) 

and 

0V  =  &  v  ,  (3.  86) 

so  that  we  can  rewrite  Eq.  's  (82)  and  (83)  as 

Ct( p,  ?)  =  k2  J  dv  CN2  fdcr  a"  8/3  J0  (r  |pv  -£v|) 


X  [  1  -  cos  (  2  p,v  )  ] 


(3.  87) 
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and 


C#  G,  ?)  =  ^“k2  1  dv  CN2  Jdc  o“*3  J0  (a|py  -?v  I  ) 

O  0 

X  [1  +  cos  (2  p,  u)  ]  .  (3.88) 

With  these  spherical  wave  propagation  theory  results  in  hand  we 
are  now  ready  to  return  to  consideration  of  adaptive  optics  systems 
performance.  At  this  point  we  are  prepared  to  undertake  the  problem  of 
evaluating  anisoplanatism  effects,  i.  e.  ,  treating  the  case  when  the  target 
direction,  0T  ,  and  the  direction  to  the  beacon,  0B  ,  are  distinct.  We  treat 
this  in  the  next  section. 


-48- 


1 . 4  System  Performance  With  Anisoplanatsim 

In  Section  1. 2  we  developed  all  of  the  l^asic  formulations  we  need 
to  study  anisoplanatism  effects  on  adaptive  optics  imaging  and  laser 
transmitter  systems,  with  and  without  random  apodization  correction. 

It  will  prove  convenient  to  gather  the  relevant  formulas  together  at 
this  point. 

For  the  adaptive  optics  imaging  system  without  random  apodization 
correction  the  relevant  formulas  are  Eq.  's  (2.  41),  (2.  42),  and  (2.  44).  We 
repeat  them  here  as 

(  I m/  o  ( f  )  =  i  (  B|  w/o  )  f  dr  (  U  iw/o  (r  +  Xf  )U|h/o  (r)  )  •  (4.  1 ) 

(R|H/n)  =  {  f  dr  |  (-U  ,n/0  (r)U|W/0  (r)  )  }  .  (4.  2) 

and 

UlHr0  (r)  =  AW  (r)  exp  {  i  k  [  R  ( 1  +  §  0T8 )  -  r  •  0T  ] 

+  i  [0(r,0T)  -0(?,0B)]+  Jt  (r,0T)  }  .  (4.3) 

For  the  adaptive  optics  imaging  system  with  random  apodization 
compensation  the  relevant  formulas  are  Eq.  's  (2.  59),  (2.  60),  and  (2.  62). 

We  repeat  them  here  as 

<  !,„(?)>  =  §<  Bin  >J  dr  <  U*IH  (7  +  xF)Uih  (r)  >  ,  (4.4) 

(B,w>  =  {  J  d 7  UtM  (r  )U,W  (7)>  }_1  (4.5) 

and 

U,  w  (  r  )  =  A  W  (  r  )  exp  (i  k  [  R  (1  +  g  0T8  )  -  r  •  0T  ] 

+  i  [0(  r,  0T  )  -  0  (  r,  0B)]  +  U(  r,  9r )  -  £(  r,  3B  )]].  (4,  6) 
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For  the  adaptive  optics  laser  transmitter  without  random 
apodization  correction  the  relevant  formulas  are  Eq.  's  (2.  68),  (2.  69),  (2.  70), 
and  (2.  72),  and  (2.  73).  We  repeat  them  here  as 

(  Glth/o  )  =  R2  (  &a»/o  )  /  (  Pltm/o  )  >  (4-  7) 

2|t 

(  ^uw/o  ^  =  H  u  m</o  uit*/o)  ,  (4.  8) 

Ulth/0  =  ^  exP  [ik  R(1+  Hr  2)  1  A  J  dr  W  (  r  ) 

X  exp  fi[0  (  r,  3T  )  -  0  (  r,  0B  )  ]  +  l  (  r ,  0T  )  ]  ,  (4.  9) 

<pi.rKi,o^  =  rdri^UiT^o(r)uLTk'o(r)i  -  (4-  10) 

and 

Ult  h/o  (  r  )  =  A  W  (  r  )  exp  f  i  k  [  O  (  r  )  +  r  •  0T  ]  -  i  0(  r,  0B)  )  .  (4.  1  1 ) 

For  the  adaptive  optics  laser  transmitter  with  random  apodization 
correction  the  relevant  formulas  are  Eq.  's  (2.  87),  (2.  88),  (2.  89),  (2.  90), 
and  (2.  105).  We  repeat  them  here  as 

(  Glt  h  )  -  R2  ( @ ltw  )  /  (  Pltw  )  >  (4-  12) 

<'PlTw)=i''uLTWulTw)  >  (4.  13) 

uLTH  =-j^exp[ikR(l+|0T8)]A  P  dr  W(r) 

X  exp  [i  [0  (  r,  0T  )  -0(  r,  0B)  ]  +  [  X(r,  0r  ) +X(r,  0B)  ]  },  (4.14) 
(PUTw>  =  J dr  |<u'LrH(r)  ULTK(r)>  ,  (4.15) 

and 

ULTW(r)  =  AW  (7)  exp  f  ik[0("r)  +  r  ■  0T  ]  -  i0(  r,  £„)  +  2(7,  0B)  }.  (4.  16) 

The  function  O  (  r  )  used  in  Eq. ' s  { 1 1  )  and  (1 6)  is  a  real  valued  function, 
defined  in  Eq.  (2.  4).  It  will  be  noted  that  Eq.  (14)  is  not  an  exact  copy  of 
Eq.  (2.  89).  Rather  it  has  been  obtained  from  Eq.  (2.  89)  by  substituting 
Eq.  's  (2.  105),  (2.  3),  and  (2.  4)  into  (2.  89)>  and  simplifying. 

With  these  formulas  in  hand  we  are  ready  to  start  a  detailed  evaluation 
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of  the  expected  performance  of  each  of  the  four  systems  under  conditions  of 
anisoplanatism.  However,  before  we  take  up  those  four  tasks  separately 
it  will  be  convenient  to  develop  some  relevant- statistical  results.  We  shall 
do  that  in  the  next  subsection. 

1.4.1  Exponential  Statistical  Functions  Associated  With  Anisoplanatism 

In  examining  the  preceding  sixteen  equations  we  can  see  that  certain 
exponential  functions  of  the  random  phase  and  log  -amplitude  will  have  to 
have  their  ensemble  average  values  determined.  These  can  be  written  as 

Eih/o  (  P  )  =  <  exp  {-  i  &S(  r  +  p,  0T)  -0(  r  +  p,  0B)  -  (  r ,  3T)  +  0(  r ,  0B)  ] 

+  [i(7+O,0T  )  +  *(?,  9t  )]  }  >  .  (4.17) 

E,M(p)  =  <  exp  f- i  [0(  r  +  p,  0t  )  -0(  r  +  p,  0B  )  -0(  r,  0T  ) +  0(  r.  0B  )  1 
+  [  i(  r  +  p,  0T  )  -  X(  r  +  p,  0B)  +  r,  0T  )  -  JL(  r,  0B)  ]  }  >  , 

(4.  18) 

Eltm/o  (  p  )  =  (  exp  f  -  i  [0(  r  +  p,  0T)  -0(  r+  p,  0B  )  -0  (  r,  0T  )  +0(  r,  0B  )  ] 

+  [X(7+pr  0T)+JKr,  0T)])>  .  (4.19) 

ELTwi  (p)  =  <  exp  { -i  [0  (  r  +  p,  0T  )  -0  (  r  +  p,  0B  )  -0(  r,  0T) +0  (  r,  0B)  ] 

+  [  *(7+p,0T)  +  .e(r  +7.  0B)  +  .i(7,eT)+.e(lr,  0b)]}>,(4.  20) 

and 

Etrwa  =  (  exP  [2X(r,Sfe)"|>  .  (4.  21) 

The  key  to  the  evaluation  of  these  expressions  lies,  of  course,  in  the 
application  of  Eq.  (2.  3  5),  and  then  in  the  use  of  Eq.  's  (3.  84),  (3.  87),  and 
(3.  88). 

Making  use  of  Eq.  (2.  35),  and  relying  on  Eq.  (3.  84)  to  allow  us  to 
drop  from  our  formulation  all  terms  involving  the  ensemble  average  of  the 
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product  of  phase  and  log  amplitude,  we  can  write 

EIM/o(p)  =  exp  { -^  <  [  0  (  r  +  p,  0T)  -0  (  r  +  p,  0B  )  -0  (  r,  0t) +0  (  r,  0B  )  ]2  > 

+  2  4  +  §<  f  LJe(r+7,eT)  -  X]+[X(7,0t)  -jt]  }2>  1  .  (4.22) 
Making  use  of  Eq.'s  (2.  37),  (3.  4),  and  (3.  5)  and  with 

?  =  0T  -  Ob  .  (4.23) 

we  can  develop  from  Eq.  (22)  the  result  that 

Eiw/o  (7)  =  exp[-2C#  (0,0)  +  2C#(7.  0)  +  2C,  (0,?)  -C,(7,  7) 

-  C,  (7,  -?)  -  cx(0,0)+  C£  (7,0)]  .  (4.24) 

Making  use  of  Eq.  's  (3.  61)  ,  (3.  82),  and  (3.  83)  we  can  rewrite  this  as 

— »  — ♦ 

E|W/o(p)  =  exp  [Slw/0  (  p  )  J  ,  (4.25) 


where 

— > 

S|w/  o  (  P  ) 


[Ci(O.O)  -cl  (7.0)]+  J^k2  j  dv 

o 


CN3  JdCT  a'*3 

o 


2  [1  +cos  (irv)]  [-1  +J0  (rrpv)] 


+  (l  +  cos[(ik^2  +  ^-S  )‘v]}[2J0((Ttfv)-Jo(e|pv-^v  I ) -J0  (® I  Pv  +7V  I)] 

2  V\ 

+  2  [1-cos  (-^-  v)  ]  [  -1  +  J0  (op, )  ]  .  (4.26) 

To  proceed  beyond  this  point  we  have  to  make  the  approximation  that  the 
turbulence  processes  are  "in  the  near  field"  so  that  we  can  approximate  the 
$  -dependent  [1  +  cos  (.  .  .)!  -term  by  2.  When  we  do  this  and  restore  the  log  - 
amplitude  covariance  dependencies  to  explicit  form,  we  get 

R  00 

Slw/o  (7)  =  [C£  (0,  0)  -Cj(7.0)]+^  ks  fdv  CN2  r d  c  o"8'3 

o  o 

X[-l+Jo  (^PvJ+Jo  (c  I  Pv  +  ^V  I)  -  4  Jo  (H  Pv  -t>»  1)3  - 

(4.  27) 
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To  evaluate  the  rr  -  integral  we  make  use  of  the  equation10 

Jdx  xwJv(ax)  =  2  u  a”  u  “ 1  T  (i  + &  v+  £tf/r  (§+  |  v  -  &  u)  , 

for  Re  (v  -1)  <  Re(u)  <5,  and  a  >  0.  (4.  28) 

Strickly  speaking  this  formula  is  not  applicable  to  our  problem  since  it 
would  involve  p  =  8/3  and  v  =  0,  for  which  the  conditions  of  applicability 
are  violated.  However,  by  arguing  that  the  -1  term  can  be  considered  to 
be  -  J0  (ax)  with  a  0,  this  formula  can  be  made  applicable  to  the  entire 
integral  of  Eq.  (27)  by  arguments  of  analytic  continuation  in  p  •  Thus  we 
obtain  the  result  that 

SlM/0  (P)  =  C  Cl  (0,0)  -  Cl  (7,0)]-  2.  905  k2  f  dv  CN8 

x[pv6/a  +  *VB/3  -  *|p,  +7V  |B/3  -I  Ipv  -?¥  |5/31-  (4-  29) 

— ♦ 

Proceeding  in  essentially  the  same  way  for  the  evaluation  of  E,  w  (  p  ) 
we  write 

E,w(p)  =  exp  { -J  (  [0  (  r  +  p,  0T)  -  0  (  r  +  7,  Sb)  -0(  r,  0T  )  +0  (  r,  0B  )]2  ) 

+i<  [*(7+0,  0T)  -X(7  +  7,9b )  +  /(7,8t)  -*(7,0b)]3>].  (4.  30) 

This  can  be  recast  in  the  form 

E|w  (p  )  =  exp  [  -2  Cf  (0,  0)  +  2C,  (7,0)  +  2C#  (0,?)  -C,(7,3)  -  C,  (7,  -7) 
+  2  Cje  (0,0)  420^(7,0)  -2CX  (0,?)  -c  x(7,?)  -Cx(p,  -*)]. 

(4.  31) 


*  The  analytic  continuation  argument  is  basically  as  follows.  Replace  the 
-8/3  power  in  Eq.  (27)  by  p  and  then  use  Eq.  (28)  to  evaluate  this  integral. 
The  result  is  an  analytic  function  of  p  which  is  clearly  valid  for  p>  -  1.  But 
since  the  quantity  in  the  square  brackets  in  Eq.  (27)  is  proportional  to  <7 2 
for  small  values  of  a,  the  integral  does  not  diverge  for  any  value  of  u 
greater  than  -3.  Accordingly  our  analytic  function  of  p  representing  the 
value  of  the  entire  integral  is  valid,  by  analytic  continuation  at  least  to 
u  =  -8/3. 
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V-  e  can  rewrite  this  as 

EIVI  (p)  =  exp  [SIVI  (p)]  ,  (4.  32) 

where 

r  8  1 6  s  R  00 

S,W(T)  =  4^  (0,0)  -C£  (0,^)]+  fdv  CN2fda  n"e/3 

0  0 

X  |2  [1+  cos(-~  v)  ]  [  -1  +  J0  (  opy )  ] 

+  f  1  +  cos  [(  rkt?3  +  )  v  ]}  [2  J0  (  rrdf )  -J0  («|  pv  -£v  | )  -J0  (o|  pv+$,  |  )3 

+  2  [1  -cos  (tt  v  )  "]  [-1  +  J0  (cP*  )  1 

+  {1  -cos  [  (ikt?2  +  |~)  v  ]  )  [2  J0  (p<Jv)-  J0(o)  Pv^v  !)  -J0(c|  Ov+«?y|)l  "jj 

=  4[CJj(0,0)-Ce(0,£)]  +  ^-k8  fdv  CNa  Jd<r  o~8/3 

x[  -1  +  Jp  (p  Pv  )+  Jo(n^v  )  'll  Jo(H  Pv -?V  ! )  '4  Jo(ct)  Pv+^v  |)3  •  (4.33) 

i  We  can  carry  out  the  ^-integration  exactly  as  we  did  before.  Doing  so  we 

obtain  the  result  that 

f  \  p 

\\  S,w(p)  =  4  [Cx(0,0)  -  Cx(0,f>)]  -  2.  905  k3  J*  dv  CNS 

i :  o 

^  x[pv6/3+^vs/3  -i!7v+^r3-i|pv  -?v  f/3]  •  (4.34) 

i'  -♦ 

j  To  evaluate  ELT  (  p)  we  note  from  a  comparison  of  Eq.  rs  (17)  and 

— ►  — ♦ 

(1 9)  that  the  expression  for  E(.TH/0  (  p  )  is  identical  to  that  for  Em/o  (  p  )  - 
Accordingly  we  can  simply  copy  Eq.  ’s  (25)  and  (29)  giving 

Eltw/o  (~o)  =  exP  [Sltm/o  ( P )  3  .  (4.35) 

*  and 
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Slt^o  (p)  =  [C^  (0,  0)  -  Cj  (  o,  0)]  -  2.  9C5  k3  f dv  CNa 

6 

x[  Pv6/a  +^vB/3  -§l7y+7,  !e/3-||py  -5v  |B/31  •  (4.  36) 

To  evaluate  ELTW1  (p)  we  start  ty  writing 

Elt  m  (  P  )  =  exp  \[-£<  [0{  r  +  p,  9T)  -0(  r  +  P,?B)  -0  (  r,  9T )  +0(  r,  9B  )]3  ) 

+4  0+4<{rA(7+p,0T )-!]+ [0(7+ p,?B )-£!  +  [£  (7, ?T )-  7] 

+  [^(7,0b)-0]]8>]|  .  (4.37) 

We  can  rewrite  this  as 

ELthi(p)  =exp[-2c,  (0,0)  +  2C#  (7,  0)  +  2C#  (0,7)  -C,  (7, 7)  -C,  (7,-7) 

-2Ci(0,0)  +  2CA(7.  0)  +  2C^(0,7)  +  C£(7,7)+Cx(7,  -7)].  (4.  38) 

Again  making  use  of  Eq.  's  (3.  82)  and  (3.  83)  we  can  rewrite  this  as 

Elt  hi  (p)  =  exp  [SLTW1  (p )  ]  ,  (4.39) 

where 

Sit  hi  (7)  =  f  dv  C*2  f  d  a  a"*3 

0  0 

X  ^2  [1+  cos  v)]  [  -1  +  J0  (  apv  )] 

+  {1+  cos  f(ik  j?2  +  %-)v  ]]  [2  J0  (  a«?v  )-J0(cr|p»'^v  |)-Jo(dPv+^t  I  )3 
+  2  [1  -cos  (^-v)]  [-1+J0  (^Pv  )  3 

+  fl  -cos[(ik02+Yv]}  [2  J0(r^)+J0(^|7v-7v|)+Jo(77v+7v|)] 

/  R  °° 

=  ^_L°ks  [ dv  C  N  2  r  do  <rfV3{[-l+  Jo(crPv)  +  Jo(o»5v)l 
0  o 

2 

-  p  cos  [(^k0s+  ^~)v]  [J0  (o|pv  -t?v | )+  J0(a|pv+tfv  I  )31  •  (4.40) 
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Here  again  we  make  the  same  approximation  as  before,  i.  e.  ,  that  we 
are  dealing  with  a  "near  field"  problem,  so  that  the  cosine -term  can  be 
approximated  by  unity.  This  allow  us  to  write 

Slt  wi  (  P  )  =  —  k2  J  dv  CN8  J  dCT  a~  0,3 

0  0 

X  [  -1  +  J0(o  Pv)  +  Jo  (ffl  Pv-«?V  \)-k  JoMpv+«?v  1)1  •  (4.  41 ) 

In  obtaining  Eq.  (29)  we  carried  out  the  same  n- integration  that  we  have  to 
perform  here.  Applying  that  result  we  can  now  write 

Slt  hi  (  o  )  =  -2.  905  ka  fdv  CN2  (p,5/3  +  ^ B/3- 1 1' pv  +  £v  p- §|pv-i>vP).  (4.42) 

o 

With  this  result  in  hand  we  now  turn  our  attention  to  the  evaluation 
of  ELtM2  ,  as  defined  in  Eq.  (21).  But  this  quantity  was  (implicitly)  evaluated 
in  going  from  Eq.  (2.  45)  to  Eq.  (2.  49).  From  that  work  we  can  see  that 

ELrK3  =  1  .  (4.43) 

With  these  results  in  hand  we  are  ready  to  take  up  the  question  of 
developing  tractable  expressions  for  system  performance  for  each  of  the 
four  types  of  systems  we  have  been  considering.  This  is  treated  in  the 
next  subsection. 

1.4.2  System  Performance  Formulas 

Before  starting  the  presentation  of  results  in  this  section  it  will  be 
convenient  to  introduce  the  following  function  which  we  will  find  is  central 
to  the  anisoplanatism  effect  for  all  cases  considered.  We  define  the  basic 
anisoplanatism  function  as 

S(pJ)  =  2.905  k2  fdv  CN3  (PyS/3+J?,6/3-i|pv+?vr/3-||pv-?v  |5/3).  (4.44) 
o 

where  we  recall  that 

p>v  =  "p  (R  -v)/R  ,  (4.45) 

$v  -  $  v  (4.  46) 
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To  evaluate  the  performance  of  an  adaptive  optics  imaging  system 
without  random  apodization  compensation  we  use  Eq.  's  (1),  (2),  (3),  (17), 
(25),  and  (29)*  It  is  easy  to  see  that  when  Eq.  (3)  is  twice  substituted  into 
Eq.  (2)  we  get 

<  B,  >  =  {  Jdr  |  A*  A  W  (?)  E,  *0  (0)  I"1 

=  {  |  A*  A  (in  Da)  EIW/0  (0)  l"1  .  (4.  47) 

From  consideration  of  Eq.  's  (25)  and  (29)  it  is  obvious  that 

E|  m/o  (0)  =  1  ,  (4.48) 

so  that 

<  EWo  >  =  {  5  A*  A  (i  tt  D8)}-1  .  (4.49) 

When  Eq.  (3)  is  twice  substituted  into  Eq.  (1)  we  get 

<2,1*,  (?)>  =  §<BIK,0>Jd?  A*  A  W(?+X?)  W  (  r  )  E,  H/0  (X? ) 

=  |  A*  A  <  B,  w/0  >  E(  h/o  (^?)  Jdr  W(?+X?)  W(?) 

=  |A*A<Blw/0)  E]Wo  (X?)  (inD3)K  (  X?)  .  (4.50) 

Now  making  use  of  Eq.  's  (2.  22),  (49),  (44),  (25),  and  (29)  we  can  write 
this  as 

<2,  h/o  (?)>  =  20l  (?)  exp  [  Cg  (0,  0)  -  Cg  (X?,  0)  -S(X?,?)]  .  (4.  51) 

It  is  insightful  to  separate  this  result  into  three  parts.  These  are:  1)  a 

— > 

diffraction  limited  term,  X0L(f ),  2)  a  random  apodization  term, 

exp  [  Cjj  (0,  0)  -  Cg  (  X f  ,  0) ] ,  which  for  most  cases  of  interest  can  be  approx  - 

imated  by  exp  [C,  (0,0)3,  and  3)  the  anisoplanatism  term, 

exp  [  -S  (  X  f  ,  d )  3  . 

For  the  evaluation  of  the  performance  of  an  adaptive  optics  imaging 
system  with  compensation  of  random  apodization  we  use  Eq.  's  (4),  (5),  (6), 
(18),  (32),  and  (34).  When  we  twice  substitute  Eq.  (6)  into  Eq.  (5)  and 
simplify  we  get 
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(4.  52) 


<  Blw  >  =  {  J*dr  4  A*A  W  (?)  E,  w  (0)  j"1 

=  {iA*A  (irrD2)  E,  M  (0)  } 

From  consideration  of  Eq.  's  (32)  and  (34)  it  is  obvious  that 

E,„  (0)  -  exp  f  4  [Cx  (0,0)  -Cj  (<),?)  ]  }  ,  (4.53) 

so  that 

<BIW  >  =  ,  §A*  A  (zttD2  )  exp  [4[CX  (0,  0)  -  Cj  •  (4.  54) 

When  Eq.  (6)  is  twice  substituted  into  Eq.  (4),  after  suitable  rearrange¬ 
ment  of  terms  we  get 

<2,w  (?)>  =  i<Blw>Jd7  A*  A  W  (?+X?)  W  (7)  EIW(X?) 

=  Ia'a  <B,m>Eiw  (Xf)  J* d?  W(?  +  X?)  W  (?) 

=  g  A*  A  <B,  w>  E(  w  ( 'T)  (JttD8)  K  (  \T)  ,  (4.55) 

where  in  writing  this  equation  we  have  suppressed  a  factor  of 
— >  ■—> 

exp  ( ik  Xf  •  9t)  as  representing  nothing  more  than  the  phase  shift  to  be 
associated  with  the  fact  that  the  object  we  are  imaging  is  displaced  an 
angular  distance  0T  from  the  z-axis.  Now  making  use  of  Eq.  (2.  22),  (54), 
(44),  (32),  and  (34)  we  can  write  this  as 

<I,M  (?)  >  =  Xou  (?)  exp  [-S  (X?,?)  ]  .  (4.56) 

In  this  case  there  is  simply  the  diffraction  limited  contribution,  XDL  (f  ), 
and  the  anisoplanalism  contribution,  exp  [-S(Xf,$  )  J  • 

To  evaluate  the  performance  of  an  adaptive  optics  laser  transmitter 
without  random  apodization  compensation  we  shall  make  use  of  Eq.  's  (7), 
(8),  (9),  (10),  (II),  (19),  (35),  and  (36).  When  we  substitute  Eq.  (II)  twice 
into  Eq.  (10)  and  suitably  simplify  we  obtain  the  result  that 

<  Plt«/o  >  =  r  A*  A  J*d?  W  (?) 

=  I  A*  A  (Jn  D2)  .  (4.  57) 
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*  i 


When  we  substitute  Eq.  (9)  twice  into  Eq.  (8),  make  a  double  integral 
of  the  product  of  integrals,  interchange  the  order  of  integration  and 
ensemble  averaging,  and  otherwise  suitably  manipulate  terms,  we  find 
that  the  result  can  be  written  as 

<^Wo>  =?A*A\-2R-2  rfdT'dr  W(r')  W  (7 )  ElTm/0  (7'- 7) 

=  iA*A\-sR-s  J]*dp  d7  W(7+T)  W(7)  Eltm/0(7) 

=  |A*A  X-2  R_s  (ir.D2)  Pd7  K(7)  Elt„/0(7)  .  (4.58) 

From  consideration  of  Eq.  's  (35)  and  (36)  we  see  that  this  can  be  rewritten  as 


<<?LTV0>=  4a*  AX-2  R-3  (|nDs)  fd7K(7)  e*p[-Cx(0,0)  +  CjU.O)-S(r.?)]- 

(4.  59) 

On  substituting  Eq.  's  (57)  and  (59)  into  Eq.  (7),  and  making  use  of  Eq.  's 
(2.  85)  and  (2.  32)  we  get  the  result  that 


<GW*,0>  =  X”  2  J d7  K  (  p )  exp  [Cx(0,  0)  -Cx(T,0)  -S(7,7)  ] 


=  t n  (  D/ >)s  K(P)  exPrq(0,0)  -Cx(p,  0)  -S(7,7)] 

Jd7  K  (p) 


=  G 


DL 


JdpK(p)  exp  [Cj  (0,  0)  -C£(7,0)  -S(p,?)] 
jd7  K  ( p)  . 


.  (4.  60) 


It  is  convenient  at  this  point  to  take  account  of  the  fact  that  for  almost  all 

cases  of  interest  the  aperture  diameter  of  the  laser  transmitter  is  very 

much  larger  than  the  range  of  the  Log -amplitude  covariance  function, 

—  ^  — ♦ 

Cg(p,0).  This  means  that  over  most  of  the  range  of  integration,  Cj(o,0) 

is  very  nearly  equal  to  zero.  This  allows  us  to  rewrite  Eq.  (60)  in  the 

approximate  form 
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(4.  61) 


<gl„/0  >  -  Gtl  «xp tc, (0, im  J^m»^p[-si^)] 

l  Jdp  K(p) 

We  can  identify  the  three  parts  of  this  result  with  1)  diffraction  limits, 

2)  random  apodization  effects,  and  3)  anisoplanatism  effects. 

Finally,  to  evaluate  the  performance  of  an  adaptive  optics  laser 
transmitter  with  compensation  for  random  apodization  we  shall  make  use 
of  Eq.  1  s  (12),  (13),  (14),  (15),  (16),  (20),  (21),  (39),  (42),  and  (43).  We 
start  by  noting  that  if  we  substitute  Eq.  (16)  into  Eq.  (15)  twice,  suitably 
manipulate  the  terms,  and  take  note  of  Eq.  (21)  we  get  the  result  that 

<PLTM  |  A*  A  f  dr  W  (7)  ELTW3  .  (4.62) 

Making  use  of  Eq.  (43)  this  reduces  to 

<  PLTH  >  =  i  A*  A  Jd7  W  (7) 

=  k  A*  A  (in  D2)  .  (4.  63) 

When  we  substitute  Eq.  (14)  twice  into  Eq.  (13),  make  a  double  integral 
of  the  product  of  integrals,  and  interchange  the  order  of  integration  and 
ensemble  averaging,  then  on  suitable  simplification  and  making  reference 
to  Eq.  (20)  we  see  that  the  result  can  be  written  as 

<^lTM  ) =  ?a*  a  x-2  r-2  JJ  d7'd7  w(7')  w(7)  e  lTW  (7' -7) 

=  |A*  a  X-2  R-2  JJ  d7  d7  w  (7  +  7)  W  (7)  elt m  (7) 

=  |A*  A  X"s  R-s  (  in  D2  )  Jd7  K  (  p  )  ELTW  (7  )  •  (4.  64) 

From  consideration  of  Eq.  ’s  (39)  and  (42)  we  see  that  this  can  be  rewritten  as 

<  0LTW  >  =  £  A*  A  X_s  R-2  (  inD2)  Jd7  K  (7)  exp  [-S<7,7)  ]  •  (4.  65) 

On  substituting  Eq.  (63)  and  (65)  into  Eq.  (12),  and  (as  before)  making  use 
of  Eq.  's  (2.  85)  and  (2.  32)  we  get  the  results  that 
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r 


<Gltm  >=  X"2  Jdp  K(  p)  exp  C-S(P,£)  ] 


£  rt(  D/X)s  .rdoK(p)exp[-S(P>«?)] 

J*  d"p  K  ( "p  ) 


_  Gb  fj  dp  K  (  p  )  exp  [  -S  (  p,  t?)  ] 

L  J  d  p  K  (  p  ) 

The  interpretation  of  this  result  is  obvious  in  the  sense  that  the  antenna 
gain  is  determined  by  two  parts,  a  diffraction  limited  part  and  an  aniso- 
planatism  part. 

1.4.3  Performance  Resume 

An  examination  of  the  results  of  the  preceding  subsection  makes  it 
clear,  as  we  might  have  expected,  that  aside  from  a  performance  factor 
reasonably  well  approximated  by  ' 

a  =  exP  Ccx  (°>°)  1  ,  (4.67) 

there  is  no  significant  difference  between  the  performance  of  adaptive  optics 
systems  designed  to  correct  for  random  apodization  and  those  not  so  designed. 
This  random  apodization  factor,  SR  h ,  is  generally  quite  close  to  unity  ,  so 
in  general  we  need  not  be  concerned  about  random  apodization. 

Ignoring  the  random  apodization  we  can  write  for  the  ensemble 
average  OTF  of  an  adaptive  optics  imaging  system 

(2  (?)  >=  X0L  exp  [  -  S  (  Xf\  &)  ]  .  (4.68) 

Similarly  ignoring  random  apodization  considerations  we  can  write  for  the 
ensemble  average  antenna  gain  of  an  adaptive  optics  laser  transmitter 

*  Implicit  in  this  is  the  assumption  that  our  propagation  path  is  such  that  the 
log -amplitude  variance,  Cjj  (0,  0),  is  less  than  one-half,  i.  e.  ,  that  scintillation 
saturation  has  not  occurred. 


(4.  66) 


-61- 


<  G  >  =  gd  j_  I  /  dp  K  (  P  )  exp  [  -  s  (  p,  0)  ] 
l  J*  d  p  K  (  p  ) 


(4.  69) 


We  recall  here  that  we  are  dealing  with  an  aperture  of  diameter  D 
and  a  wavelength  X  .  We  further  recall  that 


K(p)  = 


(fT  f  cos"1  (  p/  D  ) 
0  ,  if  p  >  D 


-  (P/D)  [  1  -  (p/D)2  ]va  if  pfiD 

(4.  70) 


SDt  (f  )  =  K  (  Xf  ) 


(4.  71) 


GdL  =  in  (D/X)2 


In  addition  we  note  that 


(4.  72) 


S<p,£)  =  2.  905  ka  Jdv  Cnb  (p*3  +  £v5/3  -  4  ft*  +  |B/3-%|'pv f/3), 

(4.  73) 


where 


pv  =  p  (R-v)/R 


(4.  74) 


1?V  =  !?  V 


(4.  75) 


These  last  nine  equations  fairly  well  sum  up  the  results  of  all  of  the 
preceding  analysis.  At  this  point  th>_  remaining  task  is  to  understand 
the  nature  of  the  S  (  p,  &)  -  function  (which  we  call  the  basic  anisoplanatism 
function)  and  its  implications  in  Eq.  's  (68)  and  (69)-  We  take  this  up  in 
the  next  section. 


1 . 5  Numerical  Results 

In  this  section  we  shall  develop  numerical  results  from  the  analytic 
expressions  presented  in  the  previous  section.  We  shall  be  interested  in 
developing  an  understanding  of  the  nature  of  the  dependence  of  anisoplanatism 
effects  on  the  various  problem  parameters,  as  evidenced  by  the  numerical 
results.  We  shall  first  undertake  an  examination  of  the  quantity  S(  p,?9), 
which  is  directly  related  to  the  anisoplanatism  effects  in  adaptive  optics 
imaging.  With  those  results  in  hand  we  will  then  take  up  consideration  of 
anisoplanatism  effect  on  antenna  gain,  (G  ),  of  an  adaptive  optics  laser  trans¬ 
mitter.  These  two  matters  are  treated  in  the  following  subsections. 

1.5.1  Adaptive  Optics  Imaging 

According  to  Eq.  (4.  68)  the  OTF  of  an  adaptive  optics  imaging  system 
as  influenced  by  anisoplanatism  effects  can  be  written  as 

<*(?>>  =  *PL  (T)  exp  [-S  (\T,  ?)]  ,  (5.1) 

— ► 

for  a  spatial  frequency  of  f  ,  an  angular  separation  between  the  direction  to 
the  beacon  and  to  the  object  to  be  imaged  of  9,  and  an  operating  wavelength  of 
X-  Here  I0  L  (f)  is  the  defraction  limited  OTF  [defined  by  Eq.  (4.71)1, 
and  the  basic  anisoplanatism  function ,  S  (Xf,$),  is  defined  in  Eq.  (4.  73). 
Replacing  the  length  Xf  by  the  length  r  we  can  write 

s(r,9)  =  2.  905  ks  J dv  CNS  f r^3[  1  -(v/R)]s/3  +  9W3  vB/3 
0 

-||r  [l-(v/R)]+£v|6/3-f  |r[l-(v/R)]  -9  v  |5/3  }  .  (5.  2) 

If  we  let  c  denote  the  cosine  of  the  angle  between  the  orientation  of  the  two 
“ ►  — ■> 

vectors,  r  and  9,  so  that 

c  =  7-  9/(  |7|  1*1  )  (5.  3) 
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then  we  can  rewrite  Eq.  (2)  as 


S(r,£)  =  2.  905  k2  Jdv  CN8  r6/3  [1  -(v/R)]E/3  +  tf6'3  vE/3 
0 

-i  {r2  [  1  -(v/R)]s  +2  rt?  [l  -(v/R)]  v  c  +  tf  2  v8  }B/6 

-i  fra[l  -(v/R)]a  -2  r  t>  [l  -(v/R)]  v  c  +  0 8  v2  J6'6  (5.4) 

It  is  convenient  at  this  point  to  consider  the  form  this  expression  takes 
for  the  two  limiting  cases  of  r  very  large  and  of  very  large. 

In  the  limiting  case  of  r  very  large  compared  to  any  value  of  0v 
(for  values  of  v  within  the  range  of  integration  and  for  which  CNS  is  not 
negligibly  small)  we  can  write 

lim  S  (7,?)  =  Sr  (r,3)  ,  (5.5) 

r  co 

where 


Sr  (r,?)  2.  905  k8  Jdv  CN8  05/s  v6'3  .  (5.  6) 

o 

The  form  of  the  integrand  in  Eq.  (6)  follows  from  Eq.  (4)  when  we  consider 
that  for  r  very  large  the  two  curly  bracket  terms  in  Eq.  (4)  each  reduce  to 
approximately  r6/3[l  -(v/R)]6^3,  and  so  together  approximately  cancel  the 
term  of  this  form  already  present  in  the  integrand.  Thus  only  the  t?5/3vB/3  - 

term  is  left  in  the  integrand  - resulting  in  Eq.  (6).  If  we  define  the  iso- 

planatic  patch  angle,  j?0  ,  by  the  equation 

0O  =  {  2.  905  ks  f  dv  CN2  v6'3  }_3'B  .  (5.7) 

o 

then  we  can  rewrite  Eq.  (6)  in  very  simple  form  as 
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Sr  (r,d)  *  («>M0)E/S  (5-  8) 

For  the  alternate  extreme  limiting  case  of  very  large  compared 
to  almost  ail  values  of  r/v  we  can  write 

lim  S  ( r,  0 )  =  Sg  ( r,  d )  ,  (5-9) 

CO 

where 


Sfl  (r,3)~  2.905  ka  Jdv  CNa  r 6/3  [  1 -(v/ R)  ]B/3  (5.10) 

0 

The  form  of  the  integrand  in  Eq.  (10)  follows  from  Eq.  (4)  when  we  consider 
that  for  #  very  large  the  two  curly  brackets  terms  in  Eq.  (4)  each  reduce, 
almost  everywhere,  to  approximately  £B/3  vs/3,  and  so  together  approximately 
cancel  the  term  of  this  form  already  present  in  the  integrand.  Thus  only 

the  r 6/3  [1  -(v/R)  ]  5/3 -term  is  left  in  the  integrand - resulting  in  Eq.  (10). 

When  we  recall  that  the  value  of  the  effective  coherence  diameter,  r0  ,  is 
given  by  the  expression11 

r0  -  {(2.  905/6.  88)  k2  /  dv  CNS  [1 -(v/R)  ]5/3}-:yB  ,  (5.11) 

0 

then  it  is  easy  to  see  that  we  can  rewrite  Eq.  (10)  in  the  rather  simple  form 

SQ  (7J)  **  6.  88  (r/r0)6/3  (5-12) 

In  conjunction  with  these  two  characteristic  propagation  parameters,  r0  and 
«?0  ,  we  introduce  the  "effective  path  length"  parameter,  s£0  ,  which  we  define 
by  the  equation 

=  r0/*0  (5-13) 
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With  these  three  parameters  in  hand  we  are  now  ready  to  numerically 

— *  — ♦ 

examine  the  behavior  of  S  (  r,  $  ). 


The  basic  anisoplanatism  function,  S  (r,i?)  is  a  function  of  two 
parameters  and  as  such  is  awkward  for  numerical  evaluation  and  for 
graphical  examination.  However,  as  can  be  seen  from  a  study  of  Eq.  (4), 
if  we  extract  a  factor  of  r  B/f6  dB/f6  from  inside  the  integrand,  the  resulting 
integral  appears  to  depend  only  on  r/t?  and  not  on  r  or  t?  separately,  i.  e.  , 
it  is  a  one  parameter  function.  This  suggests  that  we  write  in  place  of  Eq.  (4) 


S  (  r,3)  =  2.  905  ka  (  r  d)B/e  Jdv  CHs;(r/d)BM  [1  -(v/R)]6/3+  (0/r)B/s  vB/3 

0 

-i{(r/^[l-(v/R)]2+2[l-(v/R)]  v  c  +  (tf/r)  v  3  }5/s 
-i{(r/t»[l-(v/R)]S-2[l-(v/R)]  v  c  +  (tf/r)  v2  }B/S  (5.  14) 


It  is  convenient  to  make  use  of  the  propagation  parameters,  r0,  «?0  ,  and  ^ , 
to  normalize  this  expression  and  accordingly  we  write 


S(r,£)  =  2.  905  k3  (r0  0O  f^Hr/r0  )(t/#0)7*  Jdv  CNE 


0 


X  ((**<, /r)"6'6  [1  -(v/R)]6/3+(^<,/r)B/6  v5/3  ^,"6/3 
-  i  {(*?*&>  /r)’1  [1  -(v/R)  ]  8  +2  [l  -(v/r)]  v  <,-1  c  +  (#*t0/r)  vss£0“a 


1 


6/6 


[1  -(v/R)]a  2  [1  -(v/r)]  v^c+l^o/ro)^'5]6'6 

(5.  15) 


Taking  note  of  Eq.  (13)  this  can  be  rewritten  as 


S(rJ)  -  [(r  /  r0)(j?/(?0)]6/s  J  [  (r  lr0)l(dld0  )  ] 


(5.  16) 


where 
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^(Q)  =  2.  905  k2(r0  0O  i0)E/s  Jdv  CN2(f  Cf/S  [  1  -(v/R)]2+  Q“s/6vE/3^~B/a 

0 

-*{Q[l-(v/R)]2  +  2  [1  -(v/R)]  v  ai^c  +  Q'1  v2^"2  l576 

-i{Q[l-(v/R)]3  -2  [1  -{v/R^v^,-1  c  +  Q-1  v2^f3}B/6".  (5.  17) 

Further,  taking  note  of  Eq.  (13)  along  with  Eq.  (11),  this  can  be  reduced 
to  the  form 

a*{Q)  =  6.88  |  Jdv  CN2  [1-(v/R)]5/3 ]  J dv  CN2 
0  '  0 

X  li  C f/6  [1  -(v/R) ]B/3  +  O'57®  ve/a  «4>*B/3 

-4{Q[1  -(v/R)]2  +  2[1-(v/R)]  v  c  +  Q'1  v2^“3}B/6 

-  %  { Q  [  1  -(v/R)]3  -2[l-(v/R)]v£0‘l  c  +  Q”1  v3  ^f3  }B'®  •  (5.  18) 

To  see  how  this  function,  t/(Q),  behaves  we  have  developed  numerical  results 
for  the  two  cases  of  1)  a  uniform  propagation  path,  i.  e.  ,  one  in  which  the 
strength  of  turbulence,  and  thus  the  value  of  CN2,  does  not  vary  with  position 
along  the  propagation  path,  and  2)  a  ground -to -space  path  for  which  the  strength 
of  turbulence,  and  thus  CN2,  varies  with  altitude  in  the  manner  shown  in  Fig.  1. 

For  the  ground -to -space  path  it  can  be  shown  that  the  value  of  J  (Q)  is 
independent  of  the  zenith  angle,  Y.  [Of  course  for  given  values  of  r  and  i?, 
since  r0  and  0O  depend  on  the  zenith  angle,  the  appropraite  value  of  Q  will 
depend  on  Y,  but  J  (Q)  considered  as  a  function  of  Q  per  se,  does  not  depend 
on  Y.]  To  see  why  J  (Q)  is  independent  of  Y  we  first  note  that  according  to 
Eq.  's  (7)  and  (11)  t?0  is  proportional  to  [sec  (Y)]'®,B  while  r0  is  proportional 
to  [  sec  (Y)]-:yB.  Thus  it  follows  from  Eq.  (13)  that  sd0  is  proportional  to 
sec  (Y).  Examining  Eq.  (18)  we  note  first  that  all  the  (v/R) -terms  can  be 
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dropped  since  for  ground -to -space  we  can  consider  R  to  be  essentially 
infinite,  and  second  that  all  other  v -dependencies  are  actually  (v/s£o)“ 
dependencies.  Thus  it  follows  that  in  the  two  integrals  in  Eq.  (18)  [except 
for  the  two  dv -factors,  each  proportional  to  sec  (Y)  when  expressed  in 
altitude,  h,  rather  than  path  length,  v,  dependent  form,  and  thus  having 
no  composite  sec  (Y) -dependence]  the  integrand  has  no  sec  (Y)  -dependent 
terms.  This  means  that  we  can  calculate  0/  (Q)  for  the  vertical  (Y  =  0) 
ground -to -space  case  and  consider  the  results  to  apply  to  the  case  of  ground - 
to -space  propagation  at  any  zenith  angle. 

In  Fig's  2,  3,  and  3'  we  show  for  the  uniform  path  and  the  ground -to  - 

space  path,  the  form  of  J  (Q)  as  a  function  of  Q  for  the  orientation  angles 
— ►  — * 

between  r  and  £  such  that  c  takes  the  values  0.  00,  0.  25,  0.  50,  0.  75, 
and  1.  00.  In  studying  these  results  we  note  that  in  the  case  of  small  enough 
values  of  Q,  J  (Q)  behaves  approximately  as  6.  88  Cf/s  , 

lim  J  (Q)  **  6.88  Q6/6  ,  (5-19) 

Q  -*  o 

while  for  large  enough  values  of  Q,  J  (Q)  behaves  approximately  as  Q-s/6  , 

lim  J  (Q)  «  Q"s/s  (5.  20) 

Q  °° 

These  limiting  behavior  results  are  in  agreement  with  the  limiting  results 
-  ► 

for  S  ( r,  d )  presented  in  Eq.  's  (8)  and  (1  2). 

Before  passing  from  our  consideration  of  the  behavior  of  S  (  r ,  t? )  and 
of  J  (Q)  to  take  up  consideration  of  the  adaptive  optics  laser  transmitters 
antenna  gain  we  wish  to  first  look  at  approximations  for  J  (Q)  [and  thus  for 
S  ( r,  & )]  which  will  be  accurate  a  bit  closer  to  the  Q  =  1  region  than  Eq.  's 
(19)  and  (20).  To  obtain  such  results  we  start  by  noting  that  by  retaining  the 
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Q 

Figure  2.  Am  so  plana  ti  sm  Function,  e'(Q)  ,  For  a  Horizontal  Path. 

Results  are  shown  for  c  *  0.  00,  0.  25,  0.  50,  0.  75,  and  1. 00 

with  c  =  1.00  corresponding  to  the  highest  curve  and  c  =  0.00  to  the 
lowest.  The  c  =  0.  25  curve  is  so  near  coincident  with  the  c  =  0.  00 
curve  as  to  be  hardly  distinguishable. 


(D  v  w  S)  cm  <r  to 

i  i  i  ui  ♦  ♦  ♦ 

Ul  UJ  Uf  —  IlI  UJ  UJ 

Q 

Figure  3.  Ani  soplanati  sm  Function,  •/(  O ) ,  Fora  Daytime 
Vertical  Path. 

Re  suits  a  re  shown  for  c  =  0.00,  0.25,  0.  50,  0.75,  and 
1. 00  with  c  =  1. 00  corresponding  to  the  highest  curve  and  c  ' 
0.00  to  the  lower*  The  c  =0.25  curve  is  so  near  coincident 
with  the  c  =0.00  curve  as  to  be  hardly  distinguishable. 


Q 

Figure  3*  .  Ani  soplanati  sm  Function,  »*(Q),  For  a  Nighttime 
Vertical  Path. 

Results  are  shown  for  c  =  0.  00,  0.  25,  0.  50,  0.  75,  and 
1.00,  with  c  =  1.00  corresponding  to  the  highest  curve  and  r  = 
0.00  to  the  lowest.  The  c  =0.25  curve  is  so  near  coincident 
with  the  c  =  0.00  curve  as  to  be  hardly  distinguishable. 


leading  terms  in  a  power  series  expansion  we  can  write 


F  (Q)  «  Qs/S  [1  -(v/R)]B/s[l  +|  (1  -ic»)  Q"s[l  -(v/R)]“s  vs  }  , 

if  Q  »  1  ,  (5.21) 

and 

F  (Q)  ~Q'b/s  v5'3  ^o“5/3  { 1  +1-  (1  -*ca)  Q2  [l-(v/R)f  v"2^2}  , 

if  Q  «  1  ,  (5.  22) 

where  for  convenience  we  have  used  F  (Q)  defined  as  follows; 

F  (Q)  =  HQ[l-(v/R)f+2[l-(v/R)]v;£0-1  c  +  CT1  vs  £0~s  }B/e 

+  £  {Q[l  -(v/R)8  -2  [1  -(v/R)]v  =&_1  c  +  Q-1v2=&"2  }6'6.  (5.23) 

Making  use  of  Eq.  (21)  we  can  obtain  from  Eq-  (18)  the  result  that 

R  R 

J  (Q)  ~6.  88  {  J  dv  CNS  [1  -(v/R)]5'3}-1  fdv  CN3 
0  0 

X  {Q_s/s  v5'3^'3'3  4(1-1  c8)  CT7/b  [1  -(v/R)]- 1/3  V2  afo'3}, 

if  Q  >  >  1  .  (5.  24) 

Similarly,  from  Eq. 's  (22)  and  (18)  we  obtain  the  result  that 

R  R 

J  (Q)  «  6.  88  {  f  dv  CNS  [1  -(v/R)]3'3}-1  Jdv  CN2 
0  0 

X  {Qs/e[l -(v/R)]5'3 -I -  (1-i  =3)  Q7'6  [1 -(v/R)]2  v"l/3^01/3  }  , 

if  Q  «  1  .  (5.  25) 
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It  is  convenient  to  rewrite  these  results  as 


J  (Q)  ~  A  Q"s/e  [1  -  a  (  1  -  §  c2 )  Q~l/3  ]  ,  if  Q  »  1  ,  (5.  26) 


and 


J  (Q)  =  B  Cf/e  [1  -3(1  -  |  c2)  Q1/3  ]  ,  if  Q  «  1  ,  (5.  27) 


where 


A  =  6.88=£0"s/3  {  Jdv  CN2  vB/3  }  {  /  dv  CNS  [l-(v/R)]B/s  T_1  (5  28) 

0  0 

a  =  §■  =£o"l/3  (  Jdv  CNS  [1  -(v/R)]-1'3  v2  }  {  Jdv  vB/3  l"1  ,  (5.  29) 
0  o 

B  =  6.  88  ,  (5.  30) 

8  =  §•  =d01/3  f  J^v  CN2  [  1  -(v/R)]  s  v-1/3  }  {  Jdv  CN?[l-(v/R)]b/3]-i 

0  o 

(5.  31) 

Making  use  of  Eq.  's  (7),  (11),  and  (13)  we  can  rewrite  Eq-  (28)  as 

A  -  1  ,  (5.  32) 


as  we  might  have  expected  from  consideration  of  Eq.  (20). 

To  see  just  how  well  Eq.  's  (26)  and  (27)  approximate  ^  (Q)  in  the 
vicinity  of  Q  =  1 ,  in  Fig  's  4,  5,  and  5  f  we  have  replotted  the  curves  of 

Fig.  's  2,3,  and  h*  for  c  =  0.  0,  and  1.  0.  These  results  are  shown  as  the 
solid  lines.  Superimposed  on  these  we  have  plotted  J  (Q)  as  calculated  from 
Eq.  (26)  for  Q  greater  than  or  equal  to  unity  and  from  Eq.  (27)  for  Q  less 
than  or  equal  to  unity.  As  can  be  seen  these  approximations  are  quite  good 
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lE*4 


outside  the  range  from  Q  =  0.  05  to  Q  =  2.  0 . 


1.5.  2  Adaptive  Optics  .Laser  Transmitter 

With  these  results  in  hand,  particularly  the  exact  evaluations  of 
J  (Q)  presented  in  Fig.  's  2,  3,  and  3#  we  are  now  ready  to  compute  adaptive 
optics  system  performance.  The  imaging  system  OTF  results  are  so  directly 
related  to  the  of  (Q)  curves  that  there  is  no  need  for  us  to  examine  these 
results  in  detail.  Accordingly  we  turn  our  attention  to  the  laser  transmitter 
antenna  gain  results.  The  basic  result,  for  which  we  wish  to  obtain  results 
is  for  the  normalized  antenna  gain,  (  G  )/G0L,  which  we  can  write  in  accordance 
with  Eq.  (4.  69),  as 

<  G  >/Gdc  =  {  i  nD2)-1  Jdr  K  (r)  exp  [  -  S  (7,0)  ]  ,  (5.  33) 

where 

f-  [  cos"1  (r/D)  -  (r/D)  [  1  -  (r/D)a  ]1/2  }  ,  if 

K  (r)  = <  " 

LO  ,  if  r  >  D 

for  a  circular  clear  aperture  of  diameter  D.  In  obtaining  Eq.  (33)  from 
Eq.  (4.  69)  we  have  made  use  of  the  fact  that 

f  dr  K  (  r  )  =  ^ tt  D2  (5.  35) 

In  Fig.  6  we  show  results  for  the  normalized  antenna  gain  as  a  function 
of  the  normalized  aperture  diameter,  D/r0  for  a  horizontal  (CKS  constant)  path. 
Results  are  shown  for  normalized  angular  separations,  0/t%,  of  0.  25  to  1.  50  in 
steps  of  0.  25,  from  1.  5  to  4.  0  in  steps  of  0.  5,  and  values  of  5  and  6.  The  be¬ 
havior  of  these  curves  as  a  function  of  D/r0  may  be  characterized  as  a  constant 
value  of  unity  for  D/r0  much  smaller  than  unity,  then  a  simple  power -law 
decline  to  a  second  break  where  the  value  levels  off  to  a  constant  value  of 
exp  [  (  0/0o  )B/3]  . 


r  <  D 
(5.  34) 
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In  repeating  this  work  for  vertical  propagation,  daytime  and  night¬ 
time,  for  zenith  angles  from  0°  to  70°,  we  have  found  no  noticeable  devia¬ 
tion  from  the  results  presented  in  Fig.  6,  (i.  e.  ,  noticeable  when  the  figures 
are  overlayed).  Accordingly  we  suggest  that  Fig.  6  can  be  considered  to  be 
a  universal  curve. 


Figure  6.  Normalized  Laser  Transmitter  Antenna  Gain. 

Results  are  calculated  for  CH 9  constant  over  the  propagation 
path,  but  almost  exactly  match  the  results  for  day  and  for  night 
vertical  propagation  at  any  zenith  angle. 
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Appendix  to  Chapter  1 
Computer  Program  Listing 

The  programs  listed  on  the  following  pages  was  used  to  prepare 
the  numerical  data  presented  in  this  chapter. 


1 


5U 


OU 


100 


PKlJOleAtt  UPLINKS 

Kr.ftl ..  u  1  i  i &•> ; ,  y  e  i 65 ,5;  ,Zi 5 ; 

keaoku  s  ( 5  >  y  i  LRfu  *  1  mer.2  ►  r  men 3  >  1  men  i sun , E56  ► ess 
NLAL*U  X»DXyO»XX( 1U0>  »X1( 1UU)  »L0yl_U53 
—  ^«/6*UU 
E-*.i3  =  5. /3 .DU 
E/<0  ~  /  ♦  /  6  • 


AL.OUI 


AL0G(2.  ) 


1PE12 
:E  JA  = 


N  *  100 
OX  =  l./N 

sum  -  0 
3Uft2  =  0 

=  U 

.SUM*<  =  0 
DO  5U  lKi 
A  (  X  ■■ .  5  /  *D  X 

A  A  (a/  “ 

XH.1)  =  (l.~Xl**E53 
com  =  sum<Dx*xx<i; 

SUhii  -  suMiODxxxi <l; 

DUnj  =  suns*  DX*  e 1 . - x ; ~ 1 . /3 .  ■>  Kx*xi' 

bon '(  -  JOM't ;  DX*  ( 1 .  x ;  »*2*x**  t  -  1 ,  /  3 . ; 

continue: 

: ml  1 ao  =  \ z . vusxjum ; ** < -3 ,/* , >  * 
leo  =  c  2 . vo  5/ & .  ooxsurti: ; **  < -  3 . / 5 . i 
LU  -  NO/  I  MC  I  A0 

I-  053  =  LU**E53  | 

riLflTA  s*  (5./6.  1*LU*X(  -1  ./3.  ;*SUI13/Suni 
CILIA  ~  (5./0.  J*LU**(l./3,  ;*0Urt-«j/Surt2 

s-u  ei  i  e.  e  6 » 6  o ;  r  he  i  a  o »  k  u  »  l.  o  ,  alph a  »  e:e  i  a 

I  OKTIA  I  v  '  1  ML  I  AO  -  •  ,  1  NIL  1C  ,  'i  yJXy  '  NO  -  '  > 

1  .1  PE  1  / .  ‘I  •  ax  >  ■  ALP  'HA  =  '  » 1  pe  l  2 ,  <{ ,  ;;;x ,  *  i ; 

DO  200  0=1 y 165 

0  ~  a'  ♦  I  \u  05  )  /  *i .  ) 
fil'd)  ~  DL.UG ( 0 > / ALUG2 
T  LFen3  =  U**E56 
DO  150  U=1 y  5 
C  =  (K  -i  ;/*», 
son  =  o 

DU  10U  .1=1  .N 
A  “  II"  •  5  /  JKl/X 

1  ERnl  =  e  i .  -x  >  **2*0 ♦  x*x/  (  o*lo*lu  > 

IlLRliC  =  2.*<i,-x;*x*c/lu 
i erm  =  x i < i ; * r ERrtCi cxx<i>/<  r ili <n3*c u 5 o > 

1  -•  *3DU*<  Cl  ERM'  1LRM21KXCE56-*  (  lEleni  ♦  I ERM2  )  **r  56  } 

OUM  -  SUM  «•  rERM*DX 
CON  I  INI  ir- 

SlK)  =  o.oO/  £>Urt2*SUM 
y<u»IO  =  ouo 
IP  ( Ci  ♦  L  I  .  1 )  THEN 

=  .b.  00*UxxE5<!>*(  1  . -L:!;.  |  a*  ( 1 . -C*C/3.  /  *U**(  1 

LL&C 

2lK)  =  0*X(~E5*>*( 1 . -ALPHA* ( 1 .  -C*C/3. )*N**e • 1 
END  ir 
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C  PKOGKAI  UFLINKV 

ixEnL*1  3x  ,Ot  1x13/  ,LGx33)  rKf  LnhL>i)A 

I xEALkI  SLUG  x I <13 »s>  »DU<S&>»|X;iS> 

REALXO  GUI 
LOGICAL  JLDCJNll 
PI  =  3.11159265 
IFN  =  103 
CALL  OPENRDtlFN) 

00  21!  U=l»5 

RLAD t IFN )  N, (CU I) ,1=1 ,N> , (S(I, J) , 1=1, H) 

2  U  CONTINUE 

CALL  CLOSE <IFN,0> 

DO  SO  0=1, 5 
DU  30  1=1,161 

30  SLOG  \  I  *  J  )  =  ALGG(S<I*  1  »  J>/S<I,  J)  )/ALOGlG<I*-l  >/0< 

DO  10  1=1,55 

*»o  did  =  io.**(  tn~is>/&,  > 

DC ( 1 >  =  ACOS  <  0  ♦  >  — ACOS \ • 125 ) 

DU  SO  0=1,1 

SO  DC\0>  =  ACOSx  tw-1  )/S.  )-ACOS<  tu+l  >/8.  ) 

DC  ( S  ;  =  ACUS  <  .  0/  5 )  ACOS  <  1  ♦  ) 

Nl'LurS  =  12 
NPIS  =  20  0 

theta  =  o 

DO  2U0  0=1,17 
li  x  THET  A  .  LE .  1 ♦ 19  >  THEN 
THETA  =  THETA-*  ,25 
ELSE  If  K  i  HE  f'A » LL  ♦  3 . 99  ;  THEN 
THETA  =  T HE  I  A*  ♦  5 
ELSE  If  x 1  HE T A »L I ♦ 9 « 9,  T HEN 
THETA  =  THETA* 1, 

ELSE 

GO  TO  200 
END  IT 

DO  ISO  .1=1,  S3 

DXO  =  A11N1< .001, ,U01/D<1> > 

B  =  DX0**<~1 ,/NPTS) 

SU1  =  <F-X/2,J*CDX0/2.; 

X  =  0 

x'Dune  =  .false. 

lo  =  165 

DU  100  1=1, nets 

X  =  X*  DXO 
RHO  =  X*D  ( M ) 

KT56  =  x RHO*  I  HE  I A  j  *<x  x  5 . /6  •  ) 

SOX  =  SORT  ( 1-X*X; 

K  =  ATANTSOX/X) -X«SGX 
01  =  Fx-1 10/ THETA 
DU  /O  L=LO  ,  1 ,  --1 

/O  IF  UT(L; .LT.01 )  Gu  TO  00 

L.  =  1 

00  CONTINUE 

IF  <L«GT» 16*0  L=161 

LO  =  L»1 
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PROGRAM  UPLINK! 0 


i  v’Lr-iL  HI  \  16jy  f  i  \  .UfJrO  j  yl/X(  JUO  )  y  l  N  X.'  \  30  0  /  y  A  l  6  y  y  CN  (6/ 

REALMS  S<5)  y'l ERM1  ►  1  ERM2 y  i ERAS >  I  CRM r SUM t LC3 ►  CCS 
KLALXU  X  y  CJy  XX  <  300  ;  y  XI (300  >  yXO  (  300  )  y  L.O  y  L0RS3  y  LOR 
REAL  21 16C»5> 


LOGICAL*!  IDONE 

=  !Tjt/'('!yO() 

EC3  -  S./3D0 

L/V;  =  y'y/*6 


UNITE  <0yC) 

format  (•  day  or  night?*) 
Read  <  o ,  3  >  ans 

FORMAT  < Al >  • 

1"  <  AN3 . EG .  ‘  D  ‘  >  THEN 
•  IT  N  =  90 
31  N  =  10<t 

ELSE  IF  CANS.EU.  '  N*  ;  'IHEN 
EFN  =  VI 


JFN  =  IOC 
LLOb 

30  TO  ‘T 
END  IF 

CALL  OF'CNRD  t IF  N  > 


READ  t IF  N ;  N y  <  A  < I ) y 1=1 y  N  )  y  <  CNt  I >  y 1= 1  ,  N  / 
CALL  CLOSE ( IF N ► 0 > 

N  =  200 
K  =  20. EC 
dxo  =  i.E--;' 


DXUO  =  DXO 

E:  =  (  .Ul/DX0;*:x<l./N> 
31  =  1 


32  =  S 


sum  =  o 

3Uh2  =  0 
SUMS  -  0 
SUM'S  =  0 
X  =  -DXO/ 2 . 
DU  CO  1=1 yN 


X  =  X< DXO 


XOCO  =  X 
XX  (ly  =  X**ECO 

XI  <I)  =  \  1  .  “X )  **E.C3 
H  =  X*RS A<  1 ) 

DO  20  3=31 y 32 

.i.F  (A(u)  *Ll  ♦  II .  AND . A T O+l )  ♦  GE ♦  H )  30  FO  2C 

30  1 0  30 
GO  10  30 

VAL  =  LOG <  CN T  3  )  )  H_0G(FI/A\3) )  /LOG  X  A  ( 3 ' 1  )/A(3)  >  *LOG<  CN(  3+1 
C«2<l;  =  EXP < VAL) 

31  =  3 

DX1  =  DX0*L: 

Dx  <  I  >  =  LDxO  t-DXl  )/2  . 

CN2DX  =  CN2<I)KDX(I> 

sum  =  sum »  c;n2dx*xx  <  i ; 


.'CN  <  3 )  ) 


50 

60 


61 


;*o 

c 

c 

65 

66 


100 


150 

160 
20  0 


300 


sur.2  ytJfV2-» cn2.dx*xi  < i / 

SUMS  -  3UM3  .  CN2UX* <  1 .  >  W-1./J *  > 

-,un'i  =  ‘JUrt-^MJNi'UXJKci. x)w*<.*.x*~- 

DXfl  -  0X1 
CONI  1NUL 

VhCTAO^  ( 2 ♦  ?U5*SUM1*K**  1 3 . /3 .  >>**<  -3./5.  ) 
f<  0  -  <  2  ♦  y 05/6  «  UU*SUM2*R )  »*  ( "  3  •  /  5 .  > 
lo  -  no/ rut:  t  AO 
LON  -  LO/N 

cM-Hfi  “  L(  5^/6  .1  *  I  LO/R>  **<-!♦  /O  •  )  *SUM3/3UM1 

iV  lA  =  (54/6.  )*(L0/R>**<1,'"3.  )*aum/5UMj. 

WKIIE  (6.61)  THETftO»KO,LO,ALP»IA.BfclA  ^3X>|LU  =i> 

FORMAT  (  1  IIC.IAO  =‘  r  OU-..  ■**£' .  -fOA  =  ‘  ,  1PE12 .  'I///  > 

1  1PE12.*»»3X* -ALPHA  =  * , IPE12. 1 r ^LIA 

IMAX  =  X”  1 

K£5t“»”  sm?s  ARE  USED  IN  EACH  INTEGRAL. 

JSlTC  SUSTSo  H  >  »  OX  \  X  >  »  X  0  ( I  >  *K  *  A  <  1  >  *  CN2 <  X ) 

CONI XNOL 

FORMAT  <1P‘»E20.5> 

DO  200  0-1 » 165 
q  =  2 » **  (  <  0  05 )  /  T  .  > 

aiu>  ~  o 

TERMS  =  QJKixES* 

DO  150  U—1.5 
C  =  (K-l  >/*!♦ 

GUM  =  0 

rx)  iou  x-IvIMax 

v/  --  X  0  (  X  ) 

TLRM1  =  ( 1  »-X)**2*Q+X*X/  <a*LOR»I..OR> 

1  .'-X>*X*C/LOR 

1 |"  Kr/  -  X l"  ( I )  *  I EKM3  *XX  (  X  )  /  (  I  EKM3.*L0R53>  „  > 

-  ,5D0* (  ( 1 ERM1-  1 ERM2  > **E36+  <  >  tWU+  « E.KM.  > > 

3UM  —  SUii+TERM*DX<I)*GN2(I> 

CONI XNUE 

S<IO  =  6.3S3/SUM2*SUM 
Y ( J»K)  —  G(N> 

II’  (G.L-T..2)  THEN  ,.T£J,( ,  -raC/3.  )*G**(1  »/3.  >  ) 

.1.  (  J  »  K  )  »  6.33*GK5KL56Wll.  BRTA*(1. 

ELSE  £.f  .a(1  _/».  l-'HftJlYt  1  »  —  C3CC/3.  T3KC4!X«X(  "’1  ./3. ) ) 

Z(J»K)  =  U** (  -E56 )  *  (  1  *  Ad  hh* J-  » 

END  IE 

CONTINUE  , 

WRITE  (6,160)  Q,(3UO,2  U»K)»R-X^> 

I  ORMA I  <  1 1 'Ll  0 . 2 , 5  (  2X  »x_El  0  .  *.  >  ) 

CON  I XNUE 
NI-'LOTS  =  5 

“if  SpUrL™ ,  14S .  01 .  r  <  1  -  K ) .  WLOTS  > 

NI’T-.CITG  ==  o 
CONTINUE 
NF'LOTS  -  2 


1 


IHHi 


DO  'iOlJ  K-l  »S»H 
IF  (K.LXU1)  I  HEN 
NX  =  7S 
NV:  —  VO 
fclLBE 

N1  «  75 
n;:  -  vu 
END  IF 

ii  =  l^s-Nir+i 

CALL  MEL  i  FL.  ( JFN+2 »  1 65  »  01 »  T  1 1 »  K  )  r  NELO I  3  ) 
CALL  MF’LTFL\»Jr  N+/l  »N1  y  01  >  1  ►!()  »2#*NEL01i») 

CALL  MEL  ITU  UFN»-  *1  y  N2  »  0 1 1 1 1  )  t  7.t  II  »K>  »  0  ) 
NF-LO  t  S  -  U 
<IUO  CONTINUE 
END 
END 


C  PKUGKAM  UPLINK  It) 

REAL  t.U  <  160  t  '■'  (  16C/C)  »DX(300  >  ,CN2(300)  »A(6)  rCtUA) 

KEAL*8  U<C)Pl  L.PfU  »  I  EKf12 »  F ENMO »  I  £-K(i ►  SUM r E56 1 ES3 
REAL*U  X  * U  f  XX <  30 0  >  *  Xl  1 30 0  >  , X0  ( 30  0  >  »L0  »  L0K53 1  LOR 
HEAL  Z(16C»5> 

LOGICAL xl  1DONE 
EC6  —  St  / 6D 0 
ECS  =  St/  3D  U 
L/6  -  / « / 6 
WRITE  <0»5) 

C  F  ORMA  I  <  ‘  DAY  UK  NXGHIV) 

READ  (0/6;  AN3 
6  lUKflAI  (  A1 ) 

IF  (ANS.E'Ci.  '[)'  )  rULN 
IF  N  -  VO 
dF  N  =  1U‘F 

ELSE  .LI-  (  ANii.EFF.  ’  N  ’  /  1 1  IF...N 

IFN  -  91 
JFN  =  105 
ELSE 

uo  ro  *t 

END  IF 

CALL  OF'ENRD  (IF'N) 

READ  (IFN/  Nf  (  A  <  I )  t>  1=  1  /  N  >  »  (  CN  <1/  »I— 1»N) 

CALL  CLOSE  ( IF  N  r  0  ) 

N  =  2UU 
K  =  20 .EC 
DXO  =  1 . E—7 
DXOO  -•  DXO 

B  =  ( .01/DX0)**(1./N) 
dl  =  1 
dZ  -  *f 

sum  -  u 
SUM2  =  U 
suns  =  t; 

SUM'H  =  0 
X  -•  -DxO/2 . 

DU  CO  1=1  ■>  N 

X  =  X 1 DXO 
XlJ(I)  =  X 

XX  (  I  )  =  XXXECd 

XI  ( I )  =  ( 1 . "X ) **E53 
H  =  XX| \4  A  ( 1 ) 

DU  20  J=J1  /d2 

2  0  IF  (  A  ( \J  )  t  L 1  ♦  I  'I  ♦  AND  ♦  A  (  d  F- 1 )  ♦  UL  ♦  M  )  CO  I O  2C 

GO  TO  60 
GO  1C  30 

2C  DAL.  =  I..OG(CN(Jj;  *  LOG  ( H/A  ( 3)  ) /LOU  ( A  ( J  *-l  /  /  A  ( «J )  )XLOU  (  CN(  U+l )  /CN  (  d  )  ) 

CN2  < I ;  =  exp (UAL; 

30  dl  =  d 

DXl  =  DXO XB 

DXO.)  =  (DXO  *DX1  )/2. 

CNZDX  =  CN2 ( 1 ) «DX ( I ) 

SUM 1  =  SUM. I  H;N2DXXXX(I> 
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50 

60 


61 


70 

C 

C 

65 

66 


100 


ISO 

160 
20  0 


30  0 


SUM2  -  SUM2+CN20X*Xl(l> 

sun t  -  sun 3 ■ c n2dx* (  .1 .  -x >**( - -  **  : 

bU,V(  -  SUM-1*  CN2DX*  ( 1 .  ■  x  >  t  -  1  W  > 

DXO  =  0X1 
CONI  INUE 

mtTfiOlJ»  (  2  ♦ '?  0 5*SUM  1 *R** ( 8 ♦ / 3 •  >>**<  -3./S. ) 

ko  -  ( 2 . v  os/ 6 . B0*8un2*R >**(•■■  3  ♦/  j  •  > 

LO  *  KO/ THE l AO 
LOFT  “  LU/K 

L0RS3  =  LUR**ES3  ,.?uhi 

a.  |.Hc.  s  tS./6*  )*(L0/K)**(“1  ./  J»  )*oUf“Ki,  oUt  1 

UIA  =  <5«/6  «  >*(E0/K)**(1  ,/3.  )*SUM‘1/SUM,- 
WKXIE  (6.61)  lHE1A0,RU,L.0,Ai LI'HA.bt-IA 

iriAX  ^  X“  X 

rop»ll”&'  STEPS  AHC  USES  XM  EACH  INTEGRAL.  •  > 

SSl“  «^.<I».DX<I>.X.<X»««<*..CHZ<X. 
CONTINUE 

FORMAT  (1PTE20.5) 

00  200  J-l . 165 
q  -  2  ♦  **  (  ( «J  ”05  >/**♦> 

01 (J>  -  W 
TERMS  =  u**ES6 
00  150  U~1 ,5 
C  =  (K~l  )/-•• 

SUM  -  0 

DU  10 U  I-l.XMAX 

TERfll0*1  ( 1  -  "X >  **2*0 * X**/  ( a*L0K*L0R > 

f  i:  |  ■  m ■ '  us  2  ♦  *  (  1 «  X )  *X*C/  I- OK 

•rf'r.M  vi  t  I  )  *T 'EKMJ+XX ( 1 ) /  (  TERM3*E0K6«j ) 

“!.»£«£££  "wx—LSEKTEm^rERHS.-ESL, 

SUM  —  SUM+ TEK’M*OX  ( I )  *CN2 1 1 ) 

CONTINUE 

Sao  =  6  ♦  80/SUM2*SUM 

Y(J»IO  -  SOO 

1,7  !S>EW  <  1  •  H*TA*  <  1  •  -C«/3 . )  «u« 

LL<  Z(U.IO  =  Q** (  -E56)*(  1  *  -ALPHA* l 1  •  --C*C/3.  )*«**( 
END  IF 
CONTINUE 

WRITE  (6,160)  Q ,  ( S ( K )  » 2 ( U , K )  , K— 1 , >j / 

F  ORMA  I  ( 1 F  'E 1 0  <  2 , 5  (  2X  v  21:..  1 0  .  a..  )  > 

CUNT  JiNUE 

NF'LOTS  -  5 

CALL.1* MF 'E  l  FL (  JFN ,  165,01 ,  Y  (  1 , K  )  ,  NI-T..U  1  S ) 
nf-eots  =  0 

CONTINUE 
NF'LIJIS  =  2 


•LO  =' , 
2.  •>///  ) 


(1./3. ) ) 
~  i  ♦  /  3  ♦  /  y 


-86- 


DU  <*0  0  K=1  ,S»'* 

IF  (It.LU.l  >  1 HLN 
N1  =  73 
ML  ~  VO 
L.L8L 

Ml  ==  75 
N2  ~  VO 
END  IF 

II  =  :UjS  N2+1 

CALL  NPLT'F'L  <  JFN  ♦  2  r  1  d>5»  Q1 1 Y  <  1 1 K  >  ,  Nf  'LU  l  3  ) 
CALL  Mi'Ll  F  L  C  JF  N-*'l»NlyQJ.  yZ(l»IO  » 3*NPL0T8 ) 
CALL  MF'L  I  F  L  <  JI  N«<*  y  N2 1 LU  <  11  >  »  Z<  II  ,  K  )  ,  0  ) 
NF-L.01  S  =  0 
coni  j:nul 


C  PROGRAM  UF 'LINK  11 

HEALJK4  St  1  65  t  5  /  .Oil  65  /  »LU(S5)  .  K  »  L.AMDDA 
KEAL*4  SLUG  l  16b  »  5 )  ►  D  t  5b  )  »  DUl  55 )  ►  DC  (  3) 

REAL*U  bUM 
LOGICAL  .LDUNL 
PI  =  3.1415V265 

4  WRITE  (0.5) 

5  FORMAT  t*  DAY  OR  NIGHT?') 

READ  (0.6)  A  NS 

6  FORMAT  ( A1 ) 

IK  (  ANS  . CU . ' D ' )  I  MEN 
IFN  -  104 
01  N  —  11 1 

ELSE  IF  l  ANS  *  El)  ♦  '  N  '  >  I  I TEN 

IFN  =105 
JIN  =  112 
ELSE 

GO  TO  4 
END  IF 

CALL  OF-ENKD  <  IF  N  ) 

DO  2U  J=l.b 

READ l IFN)  N» (Gil) .1=1 »N) . (3(1.0) .1=1 .N) 

20  CONTINUE 

CALL  CLOSE ( IFN . 0 ) 

DO  SO  0=1.5 
DU  JO  1=1.164 

30  SLOG(I.O)  =  AL0G(3CI*l.G>/SlI.U))/AL0G(UtI«l)/l)(I) 

DO  40  n=:i»55 

4U  D(M)  =  10. **1  lM-13)/6. ) 

DC  ( 1 )  =  ACOS  (  0  .  )  -  ACOS  (  .  1 2b  ) 

DO  50  0=1.4 

JO  DC ( O )  =  riCOS (  (01 ) /•  O ♦  )  ■  - al'03 (  (0*1  )  713 .  ) 

DC  lb)  =  ACOS ( . 07 5 ) - ACOS ll.  > 

NFLOTb  =  17 
NC IS  =  200 
T  MET  A  =  0 
DO  200  0=1.17 
IF  (I  HE  I A  ♦  LL » 1  .  4?  )  I  HEN 

I HE  I A  =  T  HE  I  A  * . 25 
ELSE  IF  ( THETA » L.E .  3 .  "79 )  THEN 
THETA  =  THETA-).  5 
ELSE  IF  1 1HETA.L1 .9.9)  THEN 
( I  IE  I  A  =  T  HE  I A-)  1 . 

ELSE 

GO  10  200 
END  IF' 

DO  1 5  0  il— ;  1 . 55 

DxO  =  AfUNll .001. .001/DIN) > 

B  =  DXU*sci  -l./NI-'IS) 

son  =  o 

X  =  -DXU/2 ♦ 

1D0NE  =  .FALSE. 

LU  =  165 
DO  1U0  1=1 .NfSS 
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X  “  XU)XU 
RHU  ~  X*T.)  V  rt  / 

R  Ibo  -  %NHU*IHL IA)**<5./6. ) 

cox  -  iiuuia-x*x> 

K  =  A 1  AN  t  3UX/X )  — X*SUX 
01  =  KHU/THETA 
DU  ,0  L— L.O  » 1 » --1 

70  if  am.;  .Li  .ui  >  ou  iu  so 

L  =  1 

UU  CONTINUE 

IF  <L.GT.l69)  L—lii'H 
L0  =  L*1 
SUMl  -  0 
DU  VO  Ul-1.5 

31  -  StL.  J1  )*U3X/UtL>  )**SLOGtL»  J1 ) 

El  =  S1*R156 

II  tSl.LT.150)  3Uhl  =  SUMl  ♦LXF'(  -31  >*DC(  J1  ) 

VU  CONTINUE 

DX1  -  DX0*fc: 

1)  (  v  X»  DX1  )  *  O I  *  1 )  I HEN 

DX1  -  l.-X 
1DONL  =  .  I  RUE . 

END  II 

sun  -  sun+t H*u»sum )*<dxo+dxi >/2.*x 
IF  tlOONE)  GO  ID  110 
DXO  -  DX1 
100  CONTINUE 

WRITE,  t  0. 105)  U»M»XH)Xl 

105  FORMAT  1 ‘  INTEGRAL  NO  I  DONE  FOR  J>M>X  »*  »2I?»F10.3> 
11 0  CONTINUE 

LGtM)  -  2 . /l-'I*SUn/  C PI/ 9  •  ) 

150  CONTINUE 
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END 
END 
END 
END 
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Chapter  2 


Pulse  Laser  Backsca 
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The  analysis  begins  with  the  equation  that  describes  the  reduction 
in  the  energy  associated  with  the  laser  pulse  as  it  propagates  through  the 
scattering  medium 


E(za) 


E(zJ  exp  [-J*  dz  0(z)] 
z. 


(1) 


Here  E(z)  is  the  energy  in  the  pulse  at  altitude,  z,  and  P(z)  is  the  total 
scattering  coefficient,  which  is  assumed  to  be  a  function  of  altitude.  The 
total  energy  scattered  out  of  the  laser  beam,  Es  ,  is  then  given  by 


Es  =  E(zx)  [1  -  exp[-/  dz  e(z)]} 


(2) 


We  are  interested  in  the  case  when  E(zx)  =  1  joule.  Therefore,  only  the 
quantity,  0(z),  need  be  evaluated  to  determine  the  total  energy  scattered 
out  of  the  beam. 


We  assume  that  the  only  scattering  mechanism  is  Rayleigh  scat¬ 
tering  by  air  molecules.  Thus,  an  evaluation  of  Rayleigh  scattering  will 
lead  to  a  determination  of  S(z)  .  When  it  is  assumed  that  the  scattering 
is  incoherent,  the  angular  scattering  coefficient,  0(0),  is  given  by1 


P(0) 


^(n2-  if 

N  X4 


sin  8  0 


(3) 


where  0  is  the  scattering  angle  and  equals  90°  in  the  forward  direction, 
and  270°  is  the  backward  direction.  In  addition,  n  is  the  refractive  index, 
N  is  the  number  of  molecules  per  unit  volume,  and  X  is  the  wavelength. 
The  total  scattering  coefficient,  0,  is  obtained  by  integrating  P(0)  over  the 
solid  angle 
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0  =  JdQ  0(0) 


(4) 


TT 

We  recognize  that  J’dfl  =  2tt  J  d0  sin  0  and  Eq.  (4)  becomes 

0 


2tT3(ns-  if  r..  .  , 

nT«  Id0  8in  * 


P  = 


(5) 


The  integral  is  a  standard  form  and  equals  4/3,  therefore,  the  total 
scattering  coefficient  is  given  by  the  expression 


8n3(n2  -  l)a 

~  3N  \4 


(6) 


We  will  be  interested  in  the  scattering  mechanism  over  a  wide 
range  of  altitudes  and  for  this  reason  it  is  necessary  to  determine  how 
P  varies  with  increasing  altitude.  The  quantity  n2  -  1  is  proportional  to 
N,  the  number  of  molecules  per  unit  volume.  Therefore,  P  is  proportional 
to  N.  From  the  perfect  gas  law  N  is  proportional  to  pressure  and  inversely 
proportional  to  temperature.  Therefore,  if  P  is  known  at  one  altitude,  its 
value  at  other  altitudes  is  given  by 


0(2) 


8n3(na-  if  P(z)  T„ 
3N  \4  P0  T(z) 


(7) 


where  P0  and  T0  are  the  original  pressure  and  temperature  at  which  the 
total  scattering  coefficient  is  known,  and  P(z)  and  T(z)  are  the  pressure 
and  temperature  at  the  desired  altitude.  The  total  energy  scattered  out 
of  the  beam  is  then 


Es  =  E(zx)  \  l  -  expl  - 


8rr3(ng  -  1) 


3N  V 


,  P(z)-Q 
dz  • 


(8) 
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1 


When  it  is  assumed  that  this  energy  is  uniformly  distributed  over  a 
sphere  centered  at  an  altitude  of  (z4  +  zs)/2  the  energy  incident  on  a  receiver 
with  area  of  1  m8  is  obtained  by  multiplying  Es  by  the  solid  angle  subtended 
by  the  receiver.  The  energy  received,  ER  ,  is  then 


=  E 


'<■ 


z,  + 


(9) 


We  are  mainly  concerned  with  wavelengths  in  the  visible  regions  of 
the  spectrum  where  photon  counting  detectors  such  as  photomultipliers  are 
available.  It  is,  therefore,  more  appropriate  to  determine  the  number  of 
photons,  Np  ,  collected  by  the  receiver  aperture.  The  energy  per  photon  is 
equal  to  hc/\  where  h  is  Planck's  constant  and  c  is  the  speed  of  light.  Thus, 


Np 


Es  \ 


he 


+  z 


■7 


(10) 


Eq.  (10)  along  with  Eq.  (8)  are  the  primary  results  of  the  analysis, 
the  total  energy  scattered  out  of  the  laser  beam  and  the  total  number  of 
photons  collected  by  the  receiver  aperture  can  both  be  determined.  The 
results  for  various  wavelengths  is  the  visible  region  of  the  spectrum, 
assuming  a  one  joule  incident  pulse,  are  summarized  in  Fig.  's  1  and  2. 

To  obtain  these  results  the  following  values  of  the  parameters  were  used 


n  =  1.000293  , 

(Ha) 

N  =  2.  547  x  io2B  m-3 

* 

(Hb) 

P0  =  1013.25  mb 

(11c) 

T0  =  288.  15° K 

(lid) 
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«  10  km 
i,  =  20  km 


tx  •  15  km 
zB  =  25  km 


zl  =  20  km 
zf  =  30  km 


Fractional  Power 
Loss  into  4n 
ata  radians 

Number  Photons 
into  lm*  aperture/ 
joule  of  Incident 
Beam 

Fractional  Power 
Loss  into  4n 

ste  radians 

Number  Photons 
into  lm*  aperture/ 
joule  of  Incident 
Beam 

Fractional  Power 
Loss  into  4n 
ate  radians 

Number  Photons 
into  lm*  aperture/ 
joule  of  Incident 
Beam 

_ _ _ 

1.9642  x  10"3 

4.  3890  x  10’ 

8.  9402  x  10"‘ 

‘ 

1.  1237  x  107 

4.  0440  X  10"* 

3.2531  x  10* 

8.  1728  x  IO"3 

1.2783  x  10* 

3.  7219  x  10"* 

3.2747  x  1(T 

1. 6840  X  10'3 

9.4826  x  10* 

1.  5256  x  10'* 

2.  04  54  x  10* 

6.  9190  x  IO"3 

5.2179  x  107 

3.  1198  x  10'* 

1.  5058  x  107 

3.  1895  x  10"* 

3.  5635  X  10* 

1.4298  *  10'* 

8.9857  x  10’ 

6.4901  X  10"3 

2.6104  x  io7 

7.9667  x  10'* 

7.  1206  x  10* 

3.  5523  x  10'* 

1. 7860  x  10* 

1.  5919  *  10"* 

5.  1223  x  io7 

3.  7412  x  10"1 

1.8376  x  1CP 

1.  1032  X  IQ"1 

4.  1599  x  10* 

5.  1183  x  10"3 

1.  2352  x  10* 

Figure  l.  Atmospheric  Scattering  Results  for  Three  Different  Altitude  Ranges. 
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Figure  2.  The  Number  of  Photons  Incident  Upon  a  lm*  Receiver  Area 
for  a  One  Joule  Incident  Pulse. 

These  results  ignore  atmospheric  attenuation  between  the  receiver 
and  altitude  zx.  Vertical  propagation  is  assumed. 
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where  the  values  for  the  first  two  quantities  are  appropriate  for  air  at  sea 
level  temperature  and  pressure  conditions,  and  were  obtained  from  Ref.  1. 
In  addition  to  PQ  and  T0  ,  the  values  of  P(z)  and  T(z)  used  to  evaluate  the 
integral  in  Eq.  (8)  were  obtained  from  Ref.  2  . 

Fig.  1  summarizes  the  results  for  the  actual  wavelengths  that  w ere 
used.  Fig.  2  is  a  plot  of  Eq.  (10)  for  each  of  the  three  altitude  ranges 
considered,  indicating  that  for  a  wavelength  of  1  p  over  three  million 
photons  are  available  if  the  atmospheric  scattering  is  observed  from  20 
to  30  kilometers  in  altitude.  This  is  a  significant  number  and  can  easily 
be  detected  with  existing  devices. 
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Earl  J.  McCartney,  Optics  of  the  Atmosphere,  John  Wiley  and 
Sons  (New  York  1976),  Chapter  4. 

U.S,  Standard  Atmosphere  Supplements,  1966,  Prepared  under  the 
sponsorship  of  Environmental  Science  Services  Administration, 
National  Aeronautic s  and  Space  Administration,  and  the  United 
States  Air  Force.  Table  5.1,  Mid  -  Latitude  s,  Spring/Fall,  p  118. 
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